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Topological crystalline insulators define a new class of topological insulator phases with gapless 
surface states protected by crystalline symmetries. In this work, we present a general theory to 
classify topological crystalline insulator phases based on the representation theory of space groups. 
Our approach is to directly identify possible nontrivial surface states in a semi-infinite system with 
a specific surface, of which the symmetry property can be described by 17 two-dimensional space 
groups. We reproduce the existing results of topological crystalline insulators, such as mirror Chern 
insulators in the pm or pmm groups, Cnv topological insulators in the p4m, p31m and p6m groups, 
and topological nonsymmorphic crystalline insulators in the pg and pmg groups. Aside from these 
existing results, we also obtain the following new results: (1) there are two integer mirror Chern 
numbers (iS) in the pm group but only one (Z) in the cm or p3ml group for both the spinless and 
spinful cases; (2) for the pmm (cmm) groups, there is no topological classification in the spinless 
case but Z^ (Z^) classifications in the spinful case; (3) we show how topological crystalline insulator 
phase in the pg group is related to that in the pm group; (4) we identify topological classification 
of the p4m, p31m, and p6m for the spinful case; (5) we find topological non-symmorphic crystalline 
insulators also existing in pgg and groups, which exhibit new features compared to those in pg 
and pmg groups. We emphasize the importance of the irreducible representations for the states at 
some specific high-symmetry momenta in the classification of topological crystalline phases. Our 
theory can serve as a guide for the search of topological crystalline insulator phases in realistic 
materials. 

PACS numbers: 73.20.At, 73.43.-f, 02.20.-a 


I. INTRODUCTION 

According to the ability to conduct electric currents, 
materials can usually be classified into insulators and 
metals. Recently, researchers discovered a new class of 
materials, of which their interior cannot conduct elec¬ 
tric currents, similar to insulators, but their surfaces 
have conducting channels. Although the surface states 
can also exist in ordinary insulators, the uniqueness of 
this class of materials lies in the fact that the existence 
of surface states is related to its bulk property, which 
can be described by the mathematical theory, topology. 
Thus, these materials are dubbed “topological insulators 
(TIs)” [IHl]. The quantum Hall effect [5], discovered in 
1980, can be viewed as the first example of topologically 
non-trivial states. Under a strong magnetic field, two- 
dimensional (2D) electron gases form Landau levels and 
get localized, while at the edge, there are gapless con¬ 
ducting modes propagating along one direction, known as 
the “chiral edge state.” A recent progress in this field is 
the discovery of time-reversal (TR) invariant TIs [611T2]. 
of which the gapless edge/surface states are protected by 
TR symmetry. Different from the chiral edge states in the 
quantum Hall effect, the edge/surface states of TR invari¬ 
ant TIs contain two counter-propagating modes with op¬ 
posite spins, thus dubbed “helical edge/surface states.” 
Two branches of helical edge/surface states are related 
to each other by TR symmetry, forming the Kramers’ 
pairs. According to the Kramers’ theorem, the energy 
of one branch at the momentum k must be the same as 


that of the other branch at — k. Therefore, at the mo¬ 
menta invariant under TR, such as k = 0, two branches 
of surface states are degenerate and form a Dirac-cone¬ 
like energy dispersion. TR invariant TIs have been con¬ 
firmed in various materials, such as HgTe quantum wells 
[Hi |T3], BiSb [14], BUSes family of materials [TbHTT] . 
etc. [lang, by different experimental methods, including 
angular-resolved photon emission spectroscopy (ARPES) 
1^ . scanning tunneling microscopy (STM) 
[ 2 TH 23 ] . and transport measurements [T51I21H2S]. Since 
the topologically non-trivial edge/surface states require 
the protection from TR symmetry, these systems are also 
regarded as symmetry-protected topological states m in 
free-fermions systems. 

Based on the above discussion, we can see that degen¬ 
eracy due to symmetry plays an essential role in topo¬ 
logical states of free-fermion systems that are protected 
by symmetry. It is well known that in crystals, there 
are different types of degeneracy, including accidental 
degeneracy and essential degeneracy. The essential de¬ 
generacy arises from the presence of certain symmetries 
and is related to the representation theory of symmetry 
groups. Therefore, one may ask the following questions: 
(1) whether or not the concept of TR invariant TIs can be 
generalized to other types of symmetries, leading to new 
classes of topological states; and (2) if one can use rep¬ 
resentation theory to classify topological states that are 
protected by other types of symmetries. It was recently 
realized that crystalline symmetry can indeed protect 
new types of topological states, which are dubbed “topo- 
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logical crystalline insulators (TCIs)” [28ti46] . In partic¬ 
ular, SnTe systems are found to be a mirror-symmetry- 
protected topological state, which has recently been con¬ 
firmed by the ARPES measurement j47h^ . Various 
physical properties, including mass acquisition m , un¬ 
conventional orbital texture m, spin-filtered edge states 
[52] , quasi-particle interference [38l |53] , and interface su¬ 
perconductivity [54], have been considered. 

In this paper, we will systematically study different 
types of TCIs based on the representation theory of crys¬ 
talline symmetry groups. Instead of directly studying 
three-dimensional (3D) bulk system, we consider a 3D 
semi-infinite system with one surface and investigate non¬ 
trivial gapless surface states which cannot be gapped 
once the crystalline symmetry is preserved. Since a 3D 
semi-infinite crystal is classified by 2D space group, this 
approach allows us to classify different types of non¬ 
trivial surface states based on the representation the¬ 
ory of 2D space group. After classifying 2D topologi¬ 
cally non-trivial surface states, we can further investigate 
topological property of the 3D bulk system and identify 
the corresponding 3D bulk topological invariants. The 
details of our approach will be discussed in the following, 
but we would like first to illustrate our main results. 


A. Summary of main results 


Our main results are summarized in Table HI for 17 
2D space groups. The first column of the table is the 
group name and the second column gives the shape of 
2D surface Brillouin zone (BZ), in which “re” is for rect- 
anglar, “h” for hexagonal, “rh” for rhombic, and “s” for 
square. Different 2D surface BZs are shown in Fig.The 
next two columns list high symmetry momentum points 
(HSPs) and high symmetry momentum lines (HSLs), 
which are essential for determining topologically non¬ 
trivial surface states. Topological classifications of TCIs 
for the spinless and spinful cases are shown in the last two 
columns. Here the symbols Z and Z 2 indicate whether 
the space of ground states for the insulating phases are 
partitioned into topological sectors labeled by an integer 
or a Z 2 quantity, respectively. Z’^ (n = 2,3,4) means 
that we can define n different integer topological invari¬ 
ants. We would like to emphasize that for several space 
groups, including p4m, p31m, p6m, pgg, pmg and pAg^ 
the classification depends on which irreducible represen¬ 
tations (Irreps) of the states near the band gap belong to. 
Therefore, for these cases, we list the required condition 
of Irreps in the bracket after the classification. Below, we 
will discuss our results of TCIs for 17 2D space groups in 
several categories. 


(1) For the groups pm, p3ml and cm (Sec. Ill A), 
we have integer number of copies of topologically non¬ 
trivial surface states, which are characterized by mirror 
Chern numbers (MCNs), belonging to a Z classification, 
for both the spinless and spinful cases. Two independent 
MCNs (Z^) can be defined for the pm group, but only 


one can be defined for the p3ml and cm groups. 

(2) For the groups pmm and cmm (Sec. IIIB), MCNs 
are not allowed for the spinless case. In contrast, four 
(two) independent MCNs can be defined for the spinful 
case of pmm {cmm) group, leading to a Z^ (Z^) classifi¬ 
cation. 

(3) The pg group (Sec. IIIC) is classified by one Z 2 
topological invariant. 

(4) For the p4m, p31m and p6m (Sec. HID), the clas¬ 
sification depends on the Irreps of the states near the 
Fermi energy, and we need to illustrate the classification 
for each group separately. For the spinless case of the 
p4m group, when all the states at T and M belong to the 
doublet (or 2D) Irrep, we can define two integer topolog¬ 
ical invariants (halved mirror chirality (HMC)), giving 
rise to 1? classification. However, if any state at T or M 
belongs to one-dimensional (ID) Irreps, no topological 
invariant can be defined. For the spinful case of the p4m 
group, three topological invariants of HMC type {!?) can 
be defined, irrespective of any condition. For the p31m 
group, there are three HSPs T, K, and K' with the point 
group Csv and three HSLs with mirror symmetry. When 
the states at n HSPs (n = 2, 3) belong to doublet Irreps, 
the topological classification is Z^ with n distinct HMCs 
for both the spinful and spinless cases. If there is only 
one HSP or no HSP with the states in doublet Irreps, one 
MCN (Z) can be defined. For the p6m group, we also 
have three HSPs, in which P belongs to the group, 
K to the group, and M to the C 2 V group, and also 
three HSLs with mirror symmetry. For the spinless case, 
there are two doublet Irreps at T and one doublet Irrep 
at K. One HMC (Z classification) can be defined when 
all the states at T and K belong to doublet Irreps. The 
C 2 V group at M has four ID Irreps Ai and Bi (i = 1, 2). 
Please refer to Table. [Hjfor the definition of Ai and Bi. 
The characters of the mirror symmetry mx and my are 
the same in Irreps A^, while opposite in Irreps Bi. The 
topological classification on T-M-K depends on the Ir¬ 
reps of the bands near the Fermi surface at M. If all the 
states at M belong to Ai (or H^), another HMC can be 
defined on T-M-K. A detailed discussion about represen¬ 
tation dependent classification is given in Section. [iHEj 
For the spinful case, all the states at T and M belong to 
the doublet Irrep, while both doublet and singlet Irreps 
are possible for the states at K. Thus, three independent 
HMCs {1?) can be defined on T-K, K-M, and M-T if the 
states at K belong to doublet Irreps. Otherwise, there 
are two HMCs that can be defined on T-M and T-K-M. 

(5) For three non-symmorphic symmetry groups pmg, 
pgg, and pAg (Sec. HIE), the classification also depends 
on Irreps at the HSPs P, X, Y and M. For the spinless 
case of the pgg group, all the states at X and Y are doubly 
degenerate (doublet Irreps). The classification depends 
on ID Irreps at T and M. There are four types of different 
ID Irreps at T (M), denoted as Ai and Bi {i = 1,2). If 
all the states at T (M) belong to the Ai or Bi Irreps, one 
integer topological invariant (Z classification), dubbed 
halved glide chirality (HGC), in analogy to HMC, can be 
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FIG. 1. Schematic plot of the surface BZ of rectangular, 
square, hexagonal, and rhombic groups in (a), (b), (c), and 
(d), respectively. 


defined along the HSL X-F-Y _(X-M-Y). If both the Ai 
and Bi states exist at F and M near the band gap, one 
can no longer define HGC along these momentum lines. 
However, a Z 2 topological invariant can still be defined. 
For the spinful case, the classification is similar to that of 
the spinless case. The only difference is that now all the 
states at F and M are doubly degenerate and it depends 
on the ID Irreps at X and Y if a Z 2 or Z topological 
invariant can be defined along the momentum line F-X- 
M (F-Y-M). For the pmg group, all the states at X and 
M are doubly degenerate, and thus a HMC can always 
be defined along the line X-M (Z classification). Another 
independent Z topological invariant can be defined along 
the line X-F-Y-M if all the states at F and Y belong 
to either Ai or Bi Irreps. For the spinful case of the 
pmg group, the classification is similar and only the role 
of X, M and Y, F is switched. For the pAg group, in 
the spinless case, when the states at F and M belong to 
doublet Irreps, three independent HGCs can be defined 
on F-X, X-M and M-F. If only the states at M (F) belong 
to doublet Irrep, and all the states at F (M) belong to 
the ID Irreps Ai or Bi of point group C 4 ^, two HGCs can 
be defined on X-M (X-F) and_X-f-M (X-M-F). If both 
the Ai and Bi states exist at F (M), the HGC can only 
be defined on X-M (X-F). If the states at both F and 
M belong to singlet Irreps, no topologically non-trivial 
phase can be realized. In the spinful case, all the states 
are doubly degenerate at F and M, and thus one HGC 
can be defined on F-M. When the states at X belong to 
Ai (or Bi) of the group C 2 v^ another one HGC can be 
defined on F-X-M. 


B. Outline 

This paper is organized as follows. Section [H] is de¬ 
voted to the description of our technique approach. We 


first describe the relationship between symmetry and de¬ 
generacy and then discuss our strategy to classify topo¬ 
logically non-trivial surface states for a 3D semi-infinite 
crystal. In Sec. Uni we perform a complete analysis of 
17 2D space groups to identify possible crystal structures 
that are able to host non-trivial surface states. The topo¬ 
logical invariant is also identified for each topologically 
non-trivial case. A conclusion is drawn in Sec. lYl 


II. METHODOLOGY 
A. Symmetry and degeneracy 

In quantum mechanics, if two different eigenstates 
share the same eigenenergy, we say these two states are 
degenerate. The degeneracy is usually due to the pres¬ 
ence of symmetry, the invariance of the Hamiltonian un¬ 
der a certain operation. The relation between degeneracy 
and symmetry has been well established based on the rep¬ 
resentation theory of symmetry groups. Here we focus on 
the degeneracies induced by crystalline symmetries. Ac¬ 
cording to different types of crystalline symmetries, we 
can further classify the corresponding degeneracies into 
three types, which are discussed in details following. 

The first type of degeneracy (type-I) is due to the non¬ 
commutation relation between symmetry operations, 
which leads to the presence of high dimensional Irreps 
of the symmetry group. From the group theory, the 
eigenstates of a system can form the basis to construct 
Irreps of the symmetry group for this system and the de¬ 
generacy is equal to the dimension of the corresponding 
Irreps [ 55 ]. In an Abelian group, in which all the sym¬ 
metry operators commute with each other, only one ID 
Irreps exist. Thus, the existence of high-dimensional Ir¬ 
reps requires that at least two of the group elements do 
not commute, AB 7 ^ BA (A, B £ G, G denotes the sym¬ 
metry group of the system with the Hamiltonian H). We 
may consider a special case, of which the system has two 
symmetry operators A and B anti-commuting with each 
other {[H^A] = 0, [H^B] = 0 and {A, H} = 0). Since 
[i7. A] = 0, we can take \ 2 p) as a common eigenstate of H 
and A, i.e., B[\^) = E\^p) and A\^) = a\^). B\^) is also 
an eigenstate of the Hamiltonian with the same eigenen¬ 
ergy E since H{B\i;)) = BH\i;) = E{B\ 2 P)). At the 
same time, we have A{B\ 2 p)) = —BA\ 2 p) = —a(H|'0)), 
which indicates that Blip) and l^p) are two orthogonal 
and degenerate eigenstates of H. 

The second type of degeneracy (type-H) occurs for the 
eigenstates belonging to different Irreps of a symmetry 
group. As mentioned above, all eigenstates of a system 
can form the basis for the representations of the corre¬ 
sponding symmetry group. Let us consider two eigen¬ 
states \'ipi{k)) and \'ip 2 {k)) belonging to two different Ir¬ 
reps of the symmetry group of the Hamiltonian, where k 
is a tuning parameter and can be regarded as the momen¬ 
tum of the BZ for our case. Due to different Irreps, the 
coupling between these two states is forbidden by symme- 
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TABLE L Summary of the symmetry and topological properties of 2D space groups. SBZ: surface BZ; HSP: high symmetry 
point; HSL: high symmetry line; re: rectangular; h: hexagonal; rh: rhombic; s: square. The surface BZ of rectangular, square, 
hexagonal, and rhombic group are shown in Fig. [^(a), (b), (c), and (d), respectively. For HSL, m is for mirror symmetry and 
g for glide symmetry. Ai and Bi denote ID Irreps at the corresponding HSP while E denotes 2D Irrep. 


Group 

SBZ 

HSP 

HSL 

Topological Classihcation 

spinless 

spinful 

pm 

re 


f-X,Y-M:m^ 

= Z X Z 

z^ 

p3ml 

h 

r: C3v 

f-Mim 

z 

Z 

cm 

rh 


T-X:my 

z 

Z 

pmm 

re 

r,XM,Y-C2v 

f-X,Y-M:m^ 

r-Y,X-M:mx 

None 


cmm 

rh 

f,X,Y:C'2„ 

T-XiTTiy; 

T-Y:mx 

None 

z^ 

P9 

re 


T-X,y-M:gy 

Z 2 

Z 2 

p4m 

s 

f,M:C'4„ 

X-.Civ 

T-Ximy; 
X-M:mx] 
T-M:md 

Z2(f,M e E) 

None (f or M ^ E) 


p31m 

h 

T,K,K':C3^ 

f-K:mi 

K-K':m2 

K'-T:m3 

Z3(3 HSPs G E) 

Z2(2 HSPs G E) 

Z(general case) 

Z3(3 HSPs G E) 

I? {2 HSPs G E) 
Z(general case) 

p6m 

h 

f-.C&v 

K-.Cs^ 

M-.C2V 

f-K:mi 

f-M:m2 

K-Mims 

Z2(f G Eiii = 1, 2),K G E,U G Ai{Bi)) 
Z(f ^ Ei{i — 1, 2),K G E,M general) 
None (f 0 = 1, 2) or K 0 F) 

I?{K G E) 

Z2(K ^ E) 

P99 

re 

f,X,Y,M:C'2„ 

f-X,Y-M:s,; 
f-Y, X-M:g, 

Z2(f,MG Ai(Bi); 

TeAi{Bi),MeBi{Ai)) 

Z(f or M G Ai; f or M G Bi) 

Z 2 (general case) 

Z^{X,Y£Ai{Bi)- 

Xe Ai(Bi),Y eBi{Ai)) 
Z(X or Y G Hi;X or Y G Bi) 
Z 2 (general case) 

pmg 

re 

r,X,Y,M:C 2 „ 

f-X,Y-M:s„; 

r-Y,X-M:m^ 

Z2(Y,f G^G Y,f GSi; 

Y e Ai{Bi),f e Bi{Ai)) 

Z(general case) 

I?{%M€Ai-%M&Bi- 
X G Ai{Bi),Me Bi{Ai)) 
Z(general case) 

p4g 

s 

f,M:C'4„ 

X:C2v 

r-X:gy; 

r-M:<7d 

Z'^(f,M G E) 

Z2(M G P;, f G Ai{Bi)- P G> M) 

Z(M G S, P general; P <->■ M) 
None(f, ^AiE) 

Z\X e Ai{Bi)) 

Z(general case) 


try, indicating that the corresponding Hamiltonian must 
be block diagonal under this basis. The eigenenergies of 
these two states vary with k and at a certain k, these two 
states may cross each other. Since the coupling between 
these two states must vanish due to symmetry, this ac¬ 
cidental degeneracy at the crossing point cannot open a 
gap. 

The third type of degeneracy (type-III) originates from 
the anti-unitary symmetry operators. For a unitary sym¬ 
metry operation A, {A(l)\A 2 p) = where \ 2 p) and 10) 

are two arbitrary wave functions. In contrast, TR op¬ 
erator 0 is anti-unitary, which means that = 

((^l^))*. For spinful fermions, the TR symmetry oper¬ 
ator satisfies the relation 0^ = —1, which can lead to 
a double degeneracy in a TR invariant system accord¬ 
ing to the so-called Kramers’ theorem [55]. As discussed 
above, this degeneracy plays an essential role in protect¬ 


ing the gapless nature of surface states of TR invariant 
TIs. Moreover, magnetic systems can also possess some 
types of anti-unit ary operations, which have been dis¬ 
cussed in details in Ref. [56]. This paper mainly focuses 
on the type-I and II degeneracies, which only concern 
unitary operators. 


B. Overview of classification principles 

After understanding the above three types of degen¬ 
eracies induced by symmetry, we will next describe our 
approach for classification and show how these differ¬ 
ent types of degeneracies can lead to non-trivial surface 
states. 

Let us start from the discussion about how TR sym¬ 
metry protects the gapless nature of surface states in 
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FIG. 2. (Color online) (a) Schematic plot of a semi-infinite 
system, (b) Schematic plot of the surface states of TR in¬ 
variant TIs. r = (0, 0) and X = (tt, 0) are two TI invariant 
momenta in the surface BZ. (c) Schematic plot of the surface 
states of trivial insulators. 


TR invariant TIs and why topological surface states are 
essentially different from normal surface states, such as 
dangling bonds, when TR symmetry is preserved. Our 
discussion follows Ref. m- We consider a semi-infinite 
system with one surface along the z direction, as shown 
in Fig. [^a). The difference for the surface states be¬ 
tween TR invariant TIs and trivial insulators are shown 
in Figs.|^b) and (c). One can see that the surface states 
for TIs connect the conduction bands to valence bands 
while those for trivial insulators are isolated from either 
conduction or valence bands. This difference relies on the 
fact that all the states at r(0, 0) and X(7r, 0) in Figs.j^b) 
and (c) are doubly degenerate due to TR symmetry. For 
non-trivial surface states in Fig.[^b), the degeneracies at 
F and X point come from different surface bands (switch¬ 
ing partners) while for trivial surface states in Fig. [^c), 
degenerate states at F and X are from the same surface 
bands. Thus, the surface states in Figs, [^b) and (c) are 
essentially different. If these degeneracies are removed by 
the breaking of TR symmetry, there will be no difference 
between the surface bands in Figs. |^b) and (c). 

From this discussion, one can clearly see that the de¬ 
generacies at F and X play an essential role in protecting 
topological surface states. This physical picture can be 
naturally generalized. If one can find a similar surface 
dispersion as that in Fig.j^b) and the degeneracies at the 
F and X points are protected by other types of symmetry, 
surface states will be robust against any perturbation at 
the surface once the symmetry exists. In a crystal, degen¬ 
eracy is usually induced by crystalline symmetry. Thus, 
we will next study how different types of degeneracy due 
to crystalline symmetry can help to protect non-trivial 
surface states. 

We again consider a semi-infinite system in 3D with 
one surface as shown in Fig.j^a), which can be described 
by 2D space group instead of 3D space group since there 
is no translational symmetry along the direction normal 
to the surface. The 2D BZ is defined for this semi-infinite 
system, as shown in Fig. |^a), and can be viewed as the 
projection of 3D bulk BZ into the 2D surface. Before 
discussing surface states of this semi-infinite system, it 
is necessary to first review some basic concepts of space 
group [55] [58]. The space group G of a crystal consists of 
the symmetry operators with the form g = {r|R+T} that 


transform the crystal into itself, where r is a point group 
operation, R is a primitive translation of the correspond¬ 
ing Bravais lattice, and r is a non-primitive translation 
vector. Non-zero r appears in the compound crystals 
with more than one atoms in each primitive cell. The 
translation symmetry operators form an Abelian sub¬ 
group T. The factor group F = G/T of G with respect 
to T is isomorphic to a point group. If one can find a set 
of generators for the whole space group G, for which one 
can always choose r = 0, the group G is called a sym- 
morphic group; otherwise, it is a non-symmorphic group. 
For the subgroup T of translational operators, we can 
define the BZ in the momentum space. In the BZ, a mo¬ 
mentum k can either be preserved under an operator or 
be transformed to another one. All the symmetry opera¬ 
tions that preserve k form a subgroup of the space group 
G, usually known as a wave-vector group or a little group 
at the momentum k, denoted as Gk- The little group Gk 
also consists of translational group T and the correspond¬ 
ing factor group is denoted as Fk. The dimension of the 
Irreps of the wave-vector group Gk or the correspond¬ 
ing factor group Fk determines the degeneracy of energy 
bands at the momentum k. In Appendix we give a 
systematic review of the representation theory of space 
groups. 

By identifying the wave-vector group for each momen¬ 
tum in the BZ, we can figure out how different types of 
degeneracy occur in the whole BZ. The detailed analy¬ 
sis for 17 2D space groups will be presented in the next 
section. Here, we will first illustrate how the occurrence 
of degeneracy can lead to non-trivial surface states. We 
first look at the type-I degeneracy in a crystal. This type 
of degeneracy only occurs at some HSPs, where the corre¬ 
sponding wave-vector group consists of at least two sym¬ 
metry operations that do not commute with each other. 
To have a non-trivial surface state as that in Fig. |^b), 
we require two separate HSPs with degenerate states, 
denoted as Ki and 7^2. This degeneracy should be split 
in the line K 1 -K 2 except Ki and K 2 . The TCI due to 
type-I degeneracy, dubbed type-I TCI, can exist when 
the above conditions are satisfied. This can occur in 
both symmorphic and non-symmorphic crystals for the 
spinless case. For symmorphic crystals, since the wave 
vector group is identical to a certain point group, the 
corresponding Irreps can be directly read out from the 
character tables of the point group. By searching for 
point groups with high-dimensional Irreps, we can easily 
identify which types of crystal structures can host type- 
I TCIs. It turns out that among 17 2D space groups, 
the symmorphic groups p4m, p31m, and p6m are al¬ 
lowed for type-I TCIs. Here, we emphasize that for a 
symmorphic crystal, aside from high-dimensional Irreps, 
ID Irreps also exist. Thus, if some states belonging to 
these ID Irreps coexist with non-trivial surface states 
near the Fermi energy, the system becomes trivial. For 
non-symmorphic crystals, the situation is different and 
the eigenstates are no longer related to the representa¬ 
tions of the point group. Instead, the so-called projective 

















6 


representations [58] of point groups are required to de¬ 
scribe eigenstates in a non-symmorphic system. It turns 
out that at certain HSPs, the projective representations 
only possess high-dimensional Irreps. Correspondingly, 
all states at these momenta are at least doubly degener¬ 
ate, in contrast to the symmorphic case. Therefore, the 
type-I TCI is more robust for non-symmorphic crystals. 
The non-symmorphic groups pmg, pgg, and pAg belong 
to this case. When spin is taken into account, we need to 
consider the spinor representations (or double represen¬ 
tations) of point groups. The spinor representation can 
also be included into the theory of projective represen¬ 
tations, which is discussed in details in Appendix For 
the spinful case, the type-I TCI can exist in pmm, pmg, 
pgg, cmm, pAm^ pAg, p31m and p6m groups. 

In constrast to the type-I TCIs, the type-II degeneracy 
occurs because of multiple different Irreps. As discussed 
above, the type-II degeneracy requires a tuning parame¬ 
ter, which can be taken as the momentum k in the BZ. 
Thus, it is required that multiple Irreps should exist at 
least on a momentum line in the BZ. In the 2D space 
group, only mirror and glide symmetry can exist on a 
momentum line and allow for multiple Irreps. Thus, the 
type-II degeneracy is possible for the pm, pp, cm, pmm, 
pmp, ppp, cmm, p4m, p4p, p3ml, p31m and p6m groups. 
The TCI due to the type-II degeneracy is dubbed the 
type-II TCI. One example of the type-II TCI is the mir¬ 
ror Chern insulator [29], for which gapless surface states 
are protected by mirror symmetry. It should be empha¬ 
sized that for the groups with type-II degeneracy listed 
above, gapless surface states are not always allowed. One 
example is that mirror Chern insulators do not exist in 
the pmm for spinless fermions. We will analyze these 
different situations separately below. 

Before considering concrete 2D space groups, we sum¬ 
marize our major steps of identifying crystal space groups 
with the required symmetry for TCIs. For a given crystal, 
the wave-vector group for each momentum in the surface 
BZ of a semi-infinite system should be first identified. 
Then, based on the representation theory of wave-vector 
groups, degeneracy can be extracted for each momen¬ 
tum. If we can find two or more than two momenta with 
high-dimensional Irreps, the type-I TCI is possible. If we 
can find a line of momenta in the BZ with more than one 
Irreps, the type-II TCI can exist. It should be empha¬ 
sized that type-I and -II TCIs can co-exist in one space 
group. Our analysis of the relation between non-trivial 
surface states and degeneracy can be applied to each 2D 
space group, which will be discussed in details following. 
Although the discussion on non-trivial surface states can 
indicate topologically nontrivial phase in the bulk sys¬ 
tem based on the surface-bulk correspondence [59] . this 
approach does not directly give us topological nature of 
bulk systems, as well as bulk topological invariants. It 
also does not tell us which kinds of models can be topo¬ 
logically nontrivial. Thus, it is necessary to construct toy 
models for TCIs and compute topologically non-trivial 
surface states explicitly. We will also identify topological 



FIG. 3. (Color online) (a) Schematic plot of the ^-projection 
relation between the 3D bulk BZ and the surface BZ. The 
points in the surface BZ are marked with bar over them, (b) 
Schematic plot of the surface BZ of pm group with two my 
invariant lines high-lighted and denoted by Kyi and Ky 2 - (c) 
Schematic plot of surface states of TCI in pm group on X'- 
F-X with Col = —Cel = 1. The red and blue lines denote 
the bands in even and odd mirror subspace, respectively, (d) 
Schematic plot of the BZ of p3ml group. The MILs are de¬ 
noted by the red, blue and green lines. The three-fold invari¬ 
ant points are marked by the yellow triangles. 


invariants for each class of TCIs. 


III. CLASSIFICATION OF TCIS BASED ON 2D 
SPACE GROUPS 


Next, we will reveal our results for the classification 
of 17 2D space groups for TCI phases based on the ap¬ 
proach described above. It turns out that there are no 
TCIs in the pi, p2, p3, p4, and p6 groups. The type-II 
TCI can exist in the pm^ pg, cm, and p3ml groups for 
both single (spinless fermions) and double groups (spin¬ 
ful fermions). For the pmm and cmm groups, there is no 
TCI for the single-group case but for the double group 
case, both type-I and -II degeneracies are possible, lead¬ 
ing to a mixed type-I-II TCI phase. For the pmg, pgg, 
pAm, pAg, p31m and p6m groups, both type-I and -II 
TCIs can exist for both single and double groups. It is 
interesting to notice that type-I degeneracy is always ac¬ 
companied with type-II degeneracy for a 2D space group. 
This is because all the HSPs can be connected by mo¬ 
mentum lines that are invariant under mirror or glide 
operations. 
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A. pm, p3ml and cm groups 

We start our discussion from the systems with pm, 
p3ml and cm groups, for which only type-II degeneracy 
exists in both the spinless and spinful cases. As stated 
above, a symmetry invariant line is needed to realize the 
type-II TCI. In a 2D space group, only the mirror reflec¬ 
tion and glide symmetry operations can exist along a line. 
The most simple case is the pm group, which is a sym- 
morphic group with parallel mirror reflection axes. Let us 
assume the mirror operation is along the ^-direction, de¬ 
scribed by the notation my : {x,y) {x, —y). As shown 

in Fig. I^b), there are two mirror symmetry invariant 
momentum lines (MILs), denoted as K,yi = {kx,0) and 
i^y 2 = in the surface BZ. Along these momen¬ 

tum lines, the wave vector group is given by Cih and has 
two ID Irreps, which can be distinguished by the parity 
m± = of the mirror operator my, where / = 0 for 
the spinless case and / = 1 for the spinful case. Conse¬ 
quently, on each line, the Hamiltonian can be diagonal¬ 
ized into two blocks, which are distinguished by mirror 
parities. Thus, each block belongs to the subspace of the 
states with a definite mirror parity or m_, which is 
dubbed the mirror even or odd subspace. If two bands 
belonging to two different mirror subspaces cross each 
other along these momentum lines, the crossing point is 
protected by the mirror symmetry. 

After understanding the wave vector group of the pm 
group, we next ask what types of Fermi surfaces and en¬ 
ergy dispersions are possible for non-trivial surface states. 
Since the Fermi surface must be a closed contour on the 
surface BZ (a torus) of a bulk insulator, the Fermi surface 
always crosses the HSLs hCyi and Ky 2 with even number 
of times. There are two possibilities: (i) the Fermi sur¬ 
face intersects with each of the HSLs an even number 
of times; or (ii) it crosses both Kyi and Ky 2 with an odd 
number of times. For the case (i), we may consider an ex¬ 
ample with the Fermi surface crossing the line Kyi twice, 
as depicted by the points A and B in Fig. [^b), in which 
the solid green circle is the contour of the Fermi surface 
for surface states. The energy dispersions around these 
two points have opposite velocities for non-chiral surface 
states. We further require that the states at these two 
crossings points have opposite mirror parities so that the 
surface states can be gapless. An example of this type 
of surface energy dispersion is shown along the X'-F-X 
(Kyi) line with mirror symmetry my in Fig. [^c), where 
the red and blue lines denote surface bands with even 
and odd mirror parities, respectively. Because of oppo¬ 
site velocities, these two surface bands must cross at a 
certain momentum along the X'-f-X line. Due to op¬ 
posite mirror parities, the crossing point is protected by 
mirror symmetry and cannot open a gap. Thus, type-II 
TCI exists in the pm group. For case (ii), we may con¬ 
sider an example of the Fermi surface crossing both Kyi 
and Ky 2 once, forming an open Fermi surface, as shown 
by the dashed green line in Fig. [^b). In this case, the 
surface energy dispersion must be chiral along the lines 


Kyi and Ky 2 . This corresponds to either a layered quan¬ 
tum Hall system [60l|6T] or a Weyl semimetal [62l|63]. In 
this paper, we only focus on the insulating systems with 
non-chiral surface states in the case (i). 

Next, we discuss bulk topological invariants for the 
TCI of the pm group. As shown in Fig. I^a), we can 
always find a mirror invariant plane (MIP)7 denoted as 
Ml and M 2 in the 3D BZ, which are projected into the 
MILs Kyi and Ky 2 in the surface BZ. We can treat a MIP 
as a 2 D BZ and define the Chern number for the occupied 
states in the even (odd) subspace for a MIP as 

Ceio) = ^[ ( 1 ) 

where Fe(o)(k) is the Berry curvature of the occupied 
energy bands in the even (odd) subspace 

Fe(o)(k) = i ^ VkX«5Vk|<5). (2) 

nGocc,e(o) 

Here, \u^^) is eigenfunction of the Hamiltonian i7(k) for 
the n-th occupied energy band, and dft is an infinitesi¬ 
mal, directed area of the MIP. Since there is no coupling 
between two mirror subspaces, the Chern numbers Ce 
and Co are independent. Alternatively, one can also de¬ 
fine the mirror Chern number [64] 

CM = l{Ce-Co) (3) 

aside from the total Chern number C = \{Ce -h Co). 

The surface states for the pm group can be directly re¬ 
lated to two Chern numbers Ci and CMi (or Coi and Cei) 
defined above, where i = 1,2 denotes the MIPs Mi and 
M 2 , respectively. In the case (i), surface energy disper¬ 
sion in Fig. [^c), since there is a chiral edge mode along 
the ID line Kyi in each mirror subspace (even or odd), 
this suggests that Chern number in each mirror subspace 
is non-zero for the MIP Mi and we have Coi = — Cei = 1 
according to the bulk-edge correspondence m- This cor¬ 
responds to the non-chiral case with the total Chern num¬ 
ber Cl = C 2 = 0, but the MCN Cmi = — 1 and Cm 2 = 0. 
This type of TCI is also called “mirror Chern insulator” 
due to its non-zero MCN[29[|64|- The case (i) of the pm 
group can be characterized by two independent MCNs 
Cmi and Cm 2 [SSI, leading to the = Z x Z classifica¬ 
tion. For the case (ii), not only MCNs, but also the total 
Chern numbers in the MIPs Mi and M 2 are non-zero, 
leading to either the layered quantum Hall effect or Weyl 
semi-metals, as discussed above. 

Similar discussion of type-II TCIs can also be applied 
to the p3ml and cm groups. For the p3ml group, the 
surface BZ is shown in Fig. [^d), where bi and 62 are 
the reciprocal lattice vectors. The blue, red, and green 
lines are three equivalent MILs, and the yellow trian¬ 
gles mark the positions of three-fold-rotation-invariant 
points. Since these MILs are related to each other by 
three-fold rotations, only one independent MCN can be 
defined on the MIPs that are projected to these MILs, 


giving rise to the Z classification of the p3ml group [66] . 
The point group of f is Csv and allows for both ID and 
2D Irreps. If all the bands at T belong to the 2D Irreps, 
the crossing points between the two surface bands with 
opposite mirror parities will be pinned at the f point. 
Otherwise, the crossing points can locate at any place of 
the MILs. 

For the cm group, we consider the rhombic BZ shown 
in Fig. [^d), instead of the standard Wigner-Seitz cell. 
In this case, the mirror-symmetric line is X'-F-X for mir¬ 
ror reflection my. Different from the pm group, there is 
only one MIL in the first BZ, leading to one MCN (Z 
classification) for the cm group. 

For the pm^ p3ml and cm groups, the spinful case is 
similar to spinless case and the corresponding topological 
classifications are equivalent. 


B. pmm and cmm group 

For the pmm group, the factor group F is given by the 
point group C 2 v^ which is generated by two mirror sym¬ 
metries with perpendicular mirror reflection axes (let’s 
take them as and my). For the spinless case, the 
pmm group does not have any non-trivial surface states, 
while a mixed type-I-II TCI can exist for the spinful 
case. The essential difference between these two cases is 
the commutation relation between two mirror operators 
mx and my. Direct calculation gives mxmy = C 2 {z) and 
'^y'^x = Q{y)C 2 {z), where C 2 {z) denotes the tt rotation 
around the ^ axis and Q{y) is for the 27r rotation around 
the y axis. For the spinless case, two mirror operators 
commute with each other, [mx^my] = 0, since Q{y) = 1. 
In contrast, for the spinful case, the 27r rotation Q{y) 
gives —1, leading to {mx,my} = 0. 

Next let us check how the commutation relation be¬ 
tween mx and my makes the surface states trivial for 
the spinless case. Compared to the pm group, the pmm 
group has an additional mirror reflection symmetry. We 
consider surface bands in the line Kyi = {kx, 0) that is in¬ 
variant under my , and assume that there is a surface state 
with the momentum (/ci, 0) at the Fermi energy along the 
momentum line Kyi. The mirror operation mx requires 
that at the Fermi energy, another surface state must ex¬ 
ist with the momentum (—/ci, 0), as shown in Fig.[^ For 
spinless fermions, since mx and my commute with each 
other, these two branches of surface states at (±/ci,0) 
possess the same mirror parity under my but opposite 
velocities, as depicted in Fig. [^a). Thus, the crossing 
between these two branches at the F point (marked by 
the green circles) will open a gap, resulting the trivial 
surface states. 

In contrast, the spinful case of the pmm group is dra¬ 
matically different from the spinless case since mx and 
my anti-commute with each other. There are two conse¬ 
quences due to the anti-commutation relation. Firstly, at 
HSPs, such as F, X, Y, and M, of which the wave vector 
group includes both mx and m^, the anti-commutation 


(a) (b) 



FIG. 4. (Color online) Schematic plot of surface states of 
TCIs in the pmm group. The red and blue lines denote the 
bands with even and odd mirror parities, respectively, (a) 
The trivial phase in spinless case. The crossing marked by 
the green circle is not protected and can be gapped, (b) The 
non-trivial phase in the spinful case. 

(a) (b) 




(c) (d) 




FIG. 5. (Color online) (a) Schematic plot of non-trivial sur¬ 
face states of TR invariant TIs. (b) Schematic plot of the 
trivial surface states which are two copies of that in (a). The 
crossings marked by the green circles are not protected and 
gap can be opened, (c) Schematic plot of the non-trivial sur¬ 
face states of the pmm group in the spinful case with Cm = 1 • 
(d) Schematic plot of non-trivial surface states of the pmm 
group in the spinful case with Cm = 2. The red and blue 
lines in (c) and (d) denote the bands with even and odd mir¬ 
ror parities, respectively. 

relation yields type-I degeneracy. Secondly, the surface 
state at (/ci, 0) in the line Kyi has opposite mirror parities 
of my compared to that at (—/ci,0) for the spinful case, 
as shown in Fig. |^b), which is in sharp contrast to the 
same mirror parity for the spinless case. Since these two 
branches of surface bands in Fig. [^b) belong to different 
mirror parity subspaces, no coupling is allowed between 
these two branches. Thus, type-II degeneracy can also 
exist in this group. 

The coexistence of type-I and -II degeneracies suggests 
a Z classification of the mixed type-I-II TCI phase in the 
pmm group for the spinful case, which is different from 
the Z 2 classification of TR invariant TIs. This difference 
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can be illustrated in Fig.|^ in which we compare surface 
states for TR invariant TIs and those of TCIs in the pmm 
group. In Figs, [^a) and (b), we consider one and two 
copies of surface states for the TR invariant TIs, respec¬ 
tively. One surface state is stable while for two copies of 
surface states, the crossings (marked by the green circles) 
between them on the line connecting two TR invariant 
points [F-X in Figs.j^a) and (b)] are not protected, and 
thus a gap can be opened, leading to a trivial phase. In 
contrast, for the mixed type-I-II TCI phase of the pmm 
group, the momentum line F-X is a MIL with my sym¬ 
metry in Fig. I^c). The doubly degenerate states at the 
HSPs F and X always have opposite mirror parities of 
m^, shown by red and blue lines in Fig.|^c), respectively. 
For the case of two copies of surface states of the mixed 
type-I-II TCI phase, the crossing points along the line 
F-X are between two states with opposite mirror parities 
[Fig. [^d)] and thus protected by the type-II degeneracy. 
Therefore, in the mixed type-I-II TCI phase, any integer 
number of branches of surface bands are stable. 

All the above conclusions based on the analysis of sur¬ 
face states are consistent with the Chern number in the 
bulk systems. Detailed analysis of the total Chern num¬ 
ber and MCN in this system can be found in the Ap¬ 
pendix It turns out that the MCN for a MIP must 
be zero for the spinless case, while it can be any integer 
for the spinful case. In the surface BZ, there are four in¬ 
dependent MILs, and thus the topological classification 
of the system is = Z x Z x Z x Z for insulating sys¬ 
tems. Further analysis shows that the system possesses 
gapless Weyl points in the 3D bulk if 
number [66] . 

Just like the similarity between the pm and cm groups, 
the discussion of the cmm group is almost identical to 
that of the pmm group. The only difference is the number 
of MILs in the BZ. For the cmm group, there are only two 
MILs X'-T-X and Y'-f-Y in the rhombic BZ, as shown in 
Fig.gd). Correspondingly, the classification of the cmm 
group is 7? in the spinful case, and no topologically non¬ 
trivial phase in the spinless case. 


C. pg group 

The pg group is generated by a glide symmetry gy = 
with Tx = (1/2,0), which is a combination of 
mirror operation and nonprimitive translation. The BZ 
of the pg group is similar to that of the pm group, and 
the glide symmetry exists along the lines K,yi and Ky 2 
in Fig. I^b). The difference between the pm and pg 
group lies in the fact that the nonprimitive translation 
introduces an additional phase factor for the eigenval¬ 
ues of glide operation, compared to that of mirror op¬ 
eration. We may take the line (X'-F-X) as an ex¬ 
ample. Along this line, all the eigenstates can be char¬ 
acterized by the eigenvalues of glide operation (mirror 
operation) for the pg {pm) group and the Hamiltonian 
can be block diagonal with two blocks. For the pm 


A BA B A BA B 
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FIG. 6. (Color online) (a) Schematic plot of a lattice of the 
pm group with mirror symmetry my. The lattice constants 
in X, y, and z directions are a/2, b, and c, respectively, (b) 
Schematic plot of a lattice of the pg group, which is obtained 
by a small distortion of the lattice in (a), with glide plane 
symmetry gy. The atoms in layer A move in the y direction 
by ro, while the atoms in layer B move by —ro. The yellow 
regions in (a) and (b) denote a primitive cell of the corre¬ 
sponding lattice, (c) Schematic plot of the surface BZ of the 
lattice in (a), where Xm = (27r/a, 0) and X(^ = (—27r/a,0). 
(d) Schematic plot of the surface BZ of the lattice in (b), 
where X^ = (7r/a,0) and X^ = (—7r/a,0). 


group, these two blocks of the Hamiltonian can be la¬ 
beled by the mirror eigenvalue m± = (even or odd 
subspace). In each subspace, the block Hamiltonian is 
periodic in the BZ, so iL+(_)(X') = iL+(_)(X), where 
denotes the block Hamiltonian with the mirror 
parity + 2 -^ (“'^'^)- Thus, Chern number can be defined in 
each subspace. In contrast, for the pg group, the eigen¬ 
values of the glide operation g±{kx) = which 

is dubbed “glide parity,” depends on the momentum 
where / = 0 in the spinless case and / = 1 in the spinful 
case. When kx changes by 27r, there is an additional mi¬ 
nus sign from the phase factor As a consequence, 

g±{X') = ^zp(X) and the corresponding block Hamilto¬ 
nian in each subspace is no longer periodic. Instead, we 
have iL+(_)(X') = iL_(+)(X). Therefore, for each sub¬ 
space, the 2D BZ can not form a closed manifold to de¬ 
fine a Chern number. However, the above property of 
the glide symmetry suggests that all the bands along the 
glide symmetry invariant lines must come in pairs with 
one band of g-\-{kx) parity and the other of g-{kx) parity. 
These properties have been shown to lead to a Z 2 classi¬ 
fication of TCIs in the pg group citefang2Q15new[68] . 

Since the pg group is closely related to the pm group, 
we will discuss how TCI phases in these two groups are 
related to each other below. Here, we consider a lattice of 
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the pm group, as shown in Fig.j^a). It is a layered struc¬ 
ture stacked along the x direction, of which each layer in 
the y-z plane is a rectangular lattice, we consider the 
pm group symmetry of this lattice with mirror symme¬ 
try my : {y —y). The sites in the layers A and B are 
equivalent to each other and thus there is only one atom 
in each primitive cell for the pm group. The correspond¬ 
ing lattice constants along the x, y, z directions are a/2, 
b and c, respectively. Now let us shift the layers A and B 
in the opposite directions along the y axis with the same 
displacement ro, as shown in Fig.j^b). The mirror sym¬ 
metry is broken and changed to a glide plane symmetry 
9y — O 5 0 )}- ^ consequence, the primitive cell 

is doubled in the real space and the corresponding BZ is 
halved along the x direction, as illustrated in Figs, [^c) 
and (d). 

Next, we consider a tight-binding model with non¬ 
trivial surface states on this lattice and show how surface 
states evolve when the lattice is distorted from the pm 
group to the pg group. Let us consider the pm group case 
and assume there are four orbitals on each site, and the 
basis are denoted as {|A, 1), |A, 2), |A, 3), |A, 4)}, given 
by 


lA*) = (4) 

R 

where i = 1,2,3,4, R is the lattice vector, rA,i is the 
relative position of the atoms in one primitive cell and 
— R — ^A,i) is the corresponding atomic orbital 
wave function. Furthermore, we assume the mirror parity 
of the first two orbitals |A, 1) and |A, 2) is +1, while 
that of the other two orbitals |A, 3) and |A, 4) is — 1, 
in the spinless case. Thus, the matrix representation of 
the mirror operator in the above basis is U{my ) = Tso, 
where Tij = 0(Jj{i^j = 0,1,2,3), ai{i = 1,2,3) are 
Pauli matrices, and no is a two by two identity matrix. 
The Hamiltonian in the momentum space is given by 

Hpm(k) = [mo + to cos{kxa/2) + t'o cos{kyb) 

+to cos{kzc))]ro3 ^tsm{^{kx - 0))r3i 

sm{kzc)ro2 + t2 sm{kyb)rii, ( 5 ) 

which satisfies U(my)H{myY = iL(m^k). By 
choosing the parameters in certain regime, the Hamilto¬ 
nian can possess non-zero MCN. For example, one surface 
Dirac cone in the line X^-F-X^ is shown for the surface 
energy dispersion of this model in a slab configuration in 
Fig. mb) and two surface Dirac cones, one in XL-F-X^ 
and the other in M^-Y-M^, are shown in Fig. md). 

After distorting the lattice from the pm group symme¬ 
try to the pg group symmetry, the bases become {|A, 1), 
|A,2), |A,3), |A,4), |H,1), |H,2), |H,3), |H,4)}, and or¬ 
bital corresponding to |H, i) has the same mirror parity as 
that of |A, i). According to the matrix of mirror operator 
U{my), one can show that the glide operation exchanges 



X'g r Xg Mg Y M' X^ r Xg Mg Y M' X'g 


FIG. 7. (Color online) (a) Illustration of the BZ folding pro¬ 
cess from the lattice of pm group to pg group. The range 
of the first BZ of the lattice of pm group in the kx direction 
is [—27r/a, 27r/a], and that of pg group is [—7r/a, 7r/a]. The 
solid lines are the bands of pm group in the bulk gap, and the 
dashed lines are the bands that appear after the BZ folding. 
The band with glide even (odd) parity is drawn in red (blue) 
color. (b)The band dispersion of the model of the pm group 
in a slab configuration with one Dirac cone on X(^-f-Xm. We 
only show the surface bands in the bulk gap on one surface 
of the slab. The parameters are a = 6 = c = 1, mo = —2, 
t'o = 1.5, to — t'o — t — t' — t" — (j) — 0.4. The number of 

layers in the z direction is A = 60. (c) The band dispersion 
of the model of the pg group in a slab configuration. The 
parameters are the same with those in (b). The perturbation 
term ts = 0.2. (d)The band dispersion of the model of the 
pm group with one Dirac cone on X^^-F-Xm, and one Dirac 
cone on Mm-Y-M(,^. The parameters, which are different with 
those in (b), are mo = —0.9, to = 0.5. (e) The band disper¬ 
sion of the model of the pg group with the same parameters 
in (d) and ta = 0.2. (f) The band dispersion of the model of 
the pg group got by adding a surface potential F = 0.5 eV on 
the first layer of the slab. The other parameters are the same 
with those in (e). 
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the basis \A^i) and \B^i) as 


gy\A,i) = ±\B,i), 

gy\B,i) = ±e-^'^-'^\A,i), ( 6 ) 

where the coefficient takes + if i = 1, 2 and — if i = 3,4. 
Thus, the matrix of the glide operation is 

U'{gy) = (cos(^)cri + sin(^)(72) ® Tao. (7) 

Next we need to re-write the Hamiltonian (§ into the 
new basis of the distorted lattice. We notice that the 
nearest neighbor hopping along the x direction in the 
original lattice corresponds to the hopping between the 
adjacent A and B sublattices in the distorted lattice. Fol¬ 
lowing this rule, one can show that the Hamiltonian in 
the pg group is given by 


JJpg(k) = CTO 0 {[mo + to cos{kyb) + tg cos{kzc)]ro 3 

+ti sin(fc2c)ro2 + t2 sin(fcy6)rii} 

+ ^ [(1 + cos{kxa))ai + sm{kxa)a2] 0 Foa 




sm 4>)) 


COS 


ak^ 


cri 


+ sin (^^{kx - 0 )) sin f ^ j (72 


31 


(8) 


in the ro ^ 0 limit, where a denotes the Pauli matrices in 
the A-B sublattice space. In this limit, the Hamiltonian 
has both the glide and mirror symmetries, so the energy 
dispersion of the Hamiltonian ^ can be obtained from 
that of the Hamiltonian ([^ by BZ folding. For a non¬ 
zero ro, additional terms, such as (^Fn, can exist in 
the Hamiltonian and break mirror symmetry while 
preserving glide symmetry. 

Next, we will consider how surface states of TCIs in 
the pm group in Fig. 0 ^) evolve due to the BZ fold¬ 
ing for the pg group. This is illustrated schematically in 
Fig. B a). The surface bands between X^-X^ (X^-X^) 
will be shifted to the F-X^ (X^-F). Since this shift corre¬ 
sponds to a change of momentum kx by ^2^ ja^ the glide 
parity of the corresponding bands will be changed after 
this shift [g±{k^ + v) = 9±{kx - ^) = 
example, in Fig. |^a), the red solid line between X^-X{ 
is shifted to the blue dashed line between F-X^, where 
we have used the red and blue to represent glide parities 
g-\-{kx) and g-{kx)^ respectively. After the BZ folding, 
the surface bands have some more crossing points, as 
marked by green circles in Fig. [^a). Since two states at 
the crossing points have the same glide parity, a gap can 
be opened after we introduce additional terms to break 
the mirror symmetry down to the glide symmetry. For 
the case of one Dirac cone in Fig.j^b), we still have a sin¬ 
gle Dirac cone after BZ folding, which cannot be gapped, 
as discussed in Ref. m- However, if the original model 
© of the pm group has two Dirac cones [Fig. [^d)], two 
Dirac cones appear for the model of pg group after the 
BZ folding in Fig. [^e). But, this is a trivial phase since 


if we add an energy potential on the surface, the surface 
bands can be pushed up, and a gap can be opened for 
the whole system, as shown in Fig. [^f). Therefore, the 
topological classification of the pg group is Z 2 , which is 
consistent with results in Ref. m- 


D. p4m, pSlm, and p6m group 

The type-I and -H degeneracies coexist in the pmg, 
pgg, pAm^ pAg, p31m, p6m group in both the single and 
double groups. The essential feature of these groups 
is that the HSPs with high dimensional Irreps are al¬ 
ways connected by the mirror or glide symmetry invari¬ 
ant lines. Here we discuss three symmorphic symmetry 
groups p4m, p31m and p6m and in the next section, we 
will focus on three non-symmorphic groups pmg, pgg and 
pAg. 

For the spinless case of the groups p31m, pAm and p 6 m, 
Alexandradinata et a/. [ 66 ] have pointed out that a new 
type of topological invariant, the halved-mirror chirality, 
plays an essential role. The halved-mirror chirality is an 
integer topological invariant that characterizes the sur¬ 
face states on the half-mirror-lines (HMLs) of the surface 
BZ. The HMLs are defined as the MILs that connect two 
HSPs, at which the point groups have high dimensional 
Irreps. Correspondingly, we can also define half-mirror- 
planes (HMPs), which correspond to the planes in the 
bulk BZ that are projected onto the HMLs in the surface 
BZ. Here, we take the line F-M of the pAm group as an 
example [see the Fig. Bb) for the BZ of the pAm group]. 
The point group at the HSPs F and M is (^ 4 ^, which 
possesses one 2 D Irrep. Any state belonging to this 2 D 
Irrep is doubly degenerate at these two momenta. Mir¬ 
ror symmetry exists along the momentum line F-M, and 
thus this is a HML. The Hamiltonian on this line can be 
diagonalized into two blocks with each block labeled by 
mirror parity. Two degenerate states at F and M have 
opposite mirror parities. Because of double degeneracy 
of two states with opposite mirror parities at two ends 
of a HML, it can be proved [ 66 | that an integer topo¬ 
logical invariant, the HMC, can be defined as the differ¬ 
ence between the integral of the Berry curvature on the 
corresponding HMP in the even and odd mirror parity 
subspaces. 


X = Be - Bo, (9) 

where 

•Be(o) = T / dtdk^Fe(^o){t,k^), (10) 

Jump 

where e (o) denote the mirror even (odd) subspace and 
t denotes the momenta along the HML. This topological 
invariant determines the form of non-trivial surface states 
on the HML, e.g., F-M line. Aside from F-M, there is 
another independent HML F-X-M. Thus, for the spinless 
case, the classification of TCI in the pAm group is if 
the low-energy bands at F and M belong to the 2D Irrep 
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of C/^y group. We emphasize the importance of 2 D Irreps 
at the r and M here. If the states near the band gap at 
these two momenta belong to other ID Irreps, there is 
no topological nontrivial phase. 

For the spinful case of the pAm group, all of the spinor 
representations of the C/^y group (for F and M point) and 
C 2 v group (for X point) are 2D. Thus, all the bands at F, 
X, and M are at least doubly degenerate. Thus, we can 
identify three independent HMCs on three HMLs F-X, X- 
M, and M-F. The topological classification will change 
from 1? in the spinless case to 1? in the spinful case. We 
build a model of spinful p4m group in Appendix |C 1 [ 

The topological classification of TCIs in the p31m and 
p 6 m groups has also been discussed in Ref. [ 66 ] for the 
single-group case. One finds that the classification of 
TCI in p3Im group is 7? and in p 6 m group is Z, assum¬ 
ing that the low-energy bands at all the HSPs belong 
to the doublet Irreps of the corresponding point group. 
When the singlet Irreps are taken into consideration, the 
classifications will be different. For the p3Im group, of 
which the BZ are shown in Fig. [^c), the HSPs r, K and 
K' (belonging to group C^y) are connected by MILs. If 
only two of them belong to the 2D Irrep of group, 
such as F, and K, there are two MILs, F-K and F-K'- 
K, connecting them. Here, we notice that since C^^miV- 
K')=m(K'-K), where m(r-K') means F-K' is invariant 
under the mirror reflection operation, and the eigenvalue 
of C 3 rotation for all the ID Irreps of C^y group are I, all 
the states along the line F-K' and the line K'-K share the 
same mirror parity. Thus, the topological classification 
is Z^, characterized by two HMCs, in this case. If there 
are less than two momenta belonging to 2 D Irrep, only 
one MCN can be defined in this system, leading to a Z 
classification. For the p 6 m group, the BZ is the same 
with that of p3Im, and there are only two independent 
points F (belonging to C^y) and K (belonging to C^y) 
may have 2 D Irreps. Thus, one HMC can be defined and 
give a Z classification. The momenta F and K can be 
connected by another two independent MILs F-M and 
M-K. Whether a topological invariant with Z classifica¬ 
tion can be defined on the line F-M-K depends on the 
Irreps of the states at M. The point group of M is ( 72 ^, 
which has four ID Irreps Ai and Bi (i = 1,2). When 
all the states at M belong to Ai (or Bi), an independent 
HMC can be defined and the topological classification of 
the TCI will be Z^. We emphasize the importance of the 
ID Irreps here, and the details of this argument will be 
discussed in the Sec. |HIE] If the low-energy bands belong 
to ID Irreps at either F or K, there will be no topological 
non-trivial phase. 

For the spinful case, the classification of TCIs in the 
p31m group is the same as those of the spinless case, 
because the dimensions of the Irreps of the spinor rep¬ 
resentations of the Csy group are the same with that of 
the Irreps of the single group (spinless case). For the 
p6m group, the spinor representations at M (belonging 
to C 2 v) and F (belonging to CQy) only contain 2 D Irreps 
while those at K (belong to Csy) contain both the ID 


and 2D Irreps. Therefore, if the low-energy bands at K 
belong to the 2 D Irreps, three independent HMCs can be 
defined along the HMLs F-M, F-K, and K-M (Z^ classi¬ 
fication). When the low-energy bands at K belong to ID 
Irreps, two HMCs can be defined along the HMLs F-M 
and F-K-M (Z^ classification). TCI models of p31m and 
p6m in the spinful case are given in Appendixes |C 2 | and 
|C 3[ respectively. 


E. pmg, pgg, and pAg group 

In this section, we will consider three non-symmorphic 
groups: the pmg, pgg, and pAg groups. The main differ¬ 
ence between non-symmorphic groups and symmorphic 
groups lies in the fact that the degeneracy at certain 
HSPs in the BZ can not be described by the Irreps of 
the corresponding point groups. Instead, one needs to 
introduce the so-called projective representations, which 
will be illustrated in details following. One of the au¬ 
thors have studied the pmg group in Ref. [69|, and in 
this section, we will mainly focus on the example of the 
pgg group below. After the discussion of the pgg group, 
we will illustrate the results of the pmg and pAg groups 
for both the spinless and spinful cases. 

The generators of the pgg group include the transla¬ 
tion operators, as well as two glide-symmetry operators 
9x = {^x\Ty} and gy = with Ty = (0,1/2) and 

Tx = (1/2,0). The eigenvalues of gx are g^ = 
and of gy are g^ = where / = 0 for the 

spinless case and / = I for the spinful case. A typi¬ 
cal lattice for the pgg group is shown in Fig. [^a). In 
the BZ of the pgg group, gx symmetry exists on the mo¬ 
mentum lines F-Y and X-M, while gy symmetry appears 
on the lines F-X and Y-M. Therefore, all the states in 
these momentum lines can be labeled by the eigenvalue 
of gx or gy operator and the type-I degeneracy is possible 
along these momentum lines. For four HSPs f = (0,0), 
X = (7r,0), Y = (0,7r), and M = (7r,7r), both gx and 
gy exist, and thus the corresponding factor group F\^ of 
the wave-vector group is isomorphic to C2v^ for which 
the character table is listed in Table |n| However, the 
states at these momenta cannot be directly described 
by the Irreps of the C 2 V in Table |n| For the momen¬ 
tum k = r,X,Y,M, the representation matrix for gx{y) 
can be denoted as D^{gx{y)) = D{mx{y))^ where 

from the translational part while D{mx(y)) 
only depends on the point symmetry operation and sat¬ 
isfies the commutation relation 

D{mx)D{my) = a{gx, gy)D{my)D{mx) (II) 

with a{gx,gy) = depending 

on gx, gy and momentum k. Since [mx,my] = 0, one can 
immediately see that D{mx(y)) is not a representation 
of the point group C 2 V if tr ^ I. Instead, the so-called 
projective representations, which are discussed in details 
in Appendix!^ are required to describe the degeneracy 
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FIG. 8. (Color online)(a) Schematic plot of a lattice of the 
pgg group. The yellow region denotes a primitive cell, (b) 
The energy bands of a 2D lattice of the pgg group. Ai and 
Bi {i = 1,2) denote the Irreps of the bands at T. (c) One 
gapless surface states with all the bands belonging to Irreps 
Ai at r. (d) Two copies of gapless surface states in (c). (e 
)One gapless surface states with all the bands belonging to 
Irreps Bi at f. (f) Two copies of gapless surface states in (e). 
(g) One gapless surface states with some bands belonging to 
Irreps Ai and some bands belonging to Bi at f. (h) Two 

copies of gapless surface states in (g). The crossings marked 
by the green circles can be gapped with perturbations that 
do not break the symmetry. 


of the states. Depending on different values of a, the 
projective representation belongs to different classes. Di¬ 
rect calculation shows m = 1 — “1- This 

suggests that the states at f and M are described by the 
conventional representations of C 2 v group, denoted as Kq 
class, while those at X and Y are described by projective 
representations belonging to a non-trivial class, usually 


TABLE II. The character table of C 2 v single group. 


C2v 

E 

C 2 

rrix 

rriy 


1 

1 

1 

1 

A 2 

1 

1 

-1 

-1 

Bi 

1 

-1 

1 

-1 

B 2 

1 

-1 

-1 

1 


denoted as Ki class. For the C 2 v group, all the Irreps 
in the Kq class (conventional representations) are all ID 
(see Table 0- and thus no degeneracy occurs at F and 
M. In contrast, all the Irreps in the Ki class are 2D, 
leading to the type-II degeneracy at X and Y. The de¬ 
generacy at these two momenta can also be understood 
from the anti-commutation relation between the repre¬ 
sentation matrix D{gx) and D{gy), {D{gx), D{gy)} = 0, 
at X and Y[69]. 

The type-II degeneracy at X and Y suggests the possi¬ 
bility of type-II TCIs, which can be characterized by the 
Z 2 topological invariant [69] . Furthermore, we also notice 
that type-I degeneracy exists along the momentum lines 
F-Y, X-M, F-X and Y-M, which connect four HSPs F, X, 
Y and M. Therefore, one may ask if the mixed type-I-II 
TCIs with Z classification can exist in this system. Since 
X and Y are connected by the momenta F or M with glide 
invariant lines, we need to analyze how the glide parity 
of a state evolves along the momentum lines X-F-Y and 
X-M-Y. Let us start from the momentum X, at which all 
the states are doubly degenerate. The double degenerate 
states will be split along X-F. The two split states are 
the eigenstates of gy operators and always possess op¬ 
posite glide parities ^^and g^_^ represented by red (^^) 
and blue {g^_) in Fig. ^b), respectively. Along F-Y, all 
the states are the eigenstates of gx with the eigenvalue 
which are still labeled by the red {g^) and blue {g^) 
in Fig. [^b). It should be emphasized that the glide op¬ 
eration is gy along X-F and gx along F-Y although we 
use the same colors (red and blue) to label glide parities. 
X-F and F-Y are connected at F and it turns out that de¬ 
pending on different Irreps at F, the red (or blue) states 
along X-F can be connected to either red or blue states 
along F-Y. To see this, we need to analyze the represen¬ 
tations at f. At the F point, the glide symmetry gx or gy 
is equivalent to the mirror symmetry rrix or respec¬ 
tively. Thus, all the states are classified by the Irreps for 
the C 2 V group, as listed in the Table [TT| There are four ID 
Irreps for the C 2 V group, denoted as Ai and Bi {i = 1^2). 
The eigenstates share the same parity for rrix and my if 
they belongs to the Irreps but opposite parities for 
the Irreps Bi. Therefore, for the Irreps Ai at F, the red 
(blue) states will be connected to the red (blue) states, 
while for the Irreps Bi^ the red (blue) states will be con¬ 
nected to the blue (red) states. If all the states at the F 
point near the Fermi energy belong to the Irreps any 
integer copies of surface states are possible to exist. For 
example, we show one gapless surface states in Fig.j^c), 
of which the gapless nature is guaranteed by type-II de- 





















































14 


generacy at X and Y. In Fig.j^d), we consider two copies 
of gapless surface states and in this case, type-II degen¬ 
eracy is not enough to guarantee the gapless nature since 
they may cross at a generic momentum in the line X-F-Y. 
However, in this case, we find that all the crossings are 
always between surface bands with opposite glide pari¬ 
ties and thus cannot lead to a gap opening. Therefore, 
the topological classification in this case is Z. If all the 
states near the Fermi energy at F belong to the Irreps 
as shown in Figs. |^e) and (f) for one and two copies of 
non-trivial surface states, respectively, we also find any 
integer copies of gapless surface states are stable, indicat¬ 
ing the Z classification in this case. When the Ai states 
coexist with the Bi states at F near the Fermi energy, 
multiple copies of surface states are no longer stable. We 
show an example of one and two copies of surface states 
in Figs, [^g) and (h), respectively. One can see that for 
two copies of surface states, some crossings are between 
bands with the same glide parity, as marked by the green 
circles in Fig.j^h). A gap can be opened at these points 
and drives the system into a topologically trivial phase. 
Thus, in this case, the classification is Z 2 . Aside from the 
X-F-Y line, a similar analysis can be applied to the X- 
M-Y line and lead to a similar topological classification 
at this momentum line. The above discussion suggests 
that the classification of gapless surface states in the pgg 
group sensitively depends on the Irreps of the states at F 
and M. When the states at the F and M all belong to the 
Irreps Ai (or any integer copies of surface states are 
possible and the classification is Z x Z, in which the first 
Z is for the number of surface states along X-F-Y while 
the second is for the states along X-M-Y since these two 
momentum lines are independent. 

After identifying the classification of surface states in 
the surface BZ, we next discuss topological invariants in 
the 3D bulk systems with the pgg group. There are two 
types of topological invariants that can be defined here. 
Firstly, due to the type-II degeneracy at X and Y, a Z 2 
topological invariant can be defined, which has been dis¬ 
cussed in details for the pmg group in Ref. [69]. The dou¬ 
bly degenerate bands at X and Y can be classified into 
two groups / and II based on the eigenvalue under g^ (or 
gy). Let us consider a state on the line X-U iY-T) 
in the 3D BZ which projects onto X (or Y) [Fig. |^a)] 
satisfies iL(k)|0^) = 9xW\^) = ^+100* The 

state defined as 10^^) = where Xk is a phase 

factor, is degenerate with at X (or Y) and has the 
eigenvalue g^ = —g^ under g^. The topological invariant 
can be defined as the difference between the “doublet po¬ 
larization” Pd at X and Y, i.e., A = Pd{X) — Pd{Y) mod 
2, where = Pj - P^, and -P® = ^ / dk^{(t)^\idk, !</>£) 
{a = /, II). This is a Z 2 topological invariant, similar to 
that of TR invariant TIs. Secondly, when all the states at 
F or M belong to the ID Irreps Ai (Bi), an integer topo¬ 
logical invariant, in analogy to HMC, can be defined. Let 
us consider the line X-f-Y as an example. For a continu¬ 
ous band, if it belongs to Ai {Bi) at F, the gx glide parity 
of the band along F-Y is the same as (opposite to) the 


gy glide parity along F-X. Therefore, the bands on the 
line X-F-Y can be labeled by a single operator gy in this 
case. The bands on X-M-Y have the similar properties, 
and can be labeled by the glide parity of gy on Y-M. The 
topological invariant can be defined as 

Xi = T y dti ^ dk;,{P+ - P-), (12) 


where i = 1 labels the line X-F-Y, and i = 2 labels X-M- 
Y, ti denotes the momentum along the line i with = 0 
at X, and = I at Y, while ^2 = 0 at Y and ^2 = I at 
X. ^+(-) is the Berry curvature of the occupied bands in 
the glide parity g^(^_^^{kx) subspace of gy. This definition 
is similar to the HMC, and thus, we dub it “halved glide 
chirality” (HGC). Similarly, it can also be proved that 
when xi + X 2 is an odd number, the bulk is gapless with 
Weyl points. 

To confirm our classification, we construct a tight- 
binding model for a lattice with the pgg group sym¬ 
metry, as shown in Fig. [^a). There are four atoms in 
one primitive cell with the positions = {vx^Vy)^ yb = 
{-\-rx,ry - i), VC = {\+rx, \-ry), vd = {-rx,ry). 
The symmetry analysis below does not depend on the 
values of Vx and Vy. We assume that there are two s or¬ 
bitals on each atom for simplicity and choose the basis 
in the momentum space as 


i)(r) = —k X 


„ikR 


(13) 


where R denotes the vector of the Bravais lattice, a = 
A^B^C^D label the atoms, of which the relative position 
in a primitive cell is Tq,, i = a, 6, label the two s orbitals 
and (f denotes the Lowdin orbitals. In the above ba¬ 
sis, the Hamiltonian is periodic in the momentum space 
with the reciprocal lattice vector, i.e., i^(k) = iL(k + G), 
where G is a reciprocal lattice vector. 

Under the basis {|A, a), |H,a), |C, a), |D,a), |A,6), 
|H,6), |C, 6), |D,6)} given by (13), the explicit forms of 
the matrix U for two glide plane operations are 


U{K9x) = o-Q 

( ° 

pi(^kx kyP2 

® 0 

\ 0 


^i{kx-\-ky) l 2 Q 

0 0 

0 0 

Q ^'i{~kx—ky)/2 


and 


U{K9y) = CTo 

( ° 

^ d(k.-ky)l2 

\ 0 


Q gd kx ky)j2 

0 0 

0 0 

^i{kx-\-ky) I2 Q 


0 \ 

0 

gi( —/Ccc+Zcy )/2 

0 J 

(14) 


0 \ 

gi( kx ~\~ky) / 2 

0 

0 / 

(15) 
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where (Tq is a two-by-two identity matrix. where 

The Hamiltonian in momentum space, which satis¬ 
fies the constraints enforced by symmetries = 

U{'k,g^(^y))H{k)U{k,g^(^y-jy, is 


H = 


Ha Hab\ 

hL Hb J ’ 


(16) 


Ha(k) 



0 


/2a(e*^^“+e ^kx^ 


\ 


+ e*^i»e-*'=»') 

0 


/2a(e-*^=“ +e*^=“e-*''-) 


ha{e~^^^ 

a _j_ g*02a ^ikx ^ 

/3a + /Le'"^ 

0 



1 

B 




0 

/ 


+ [-l + 8(3 + cos{kx) - cos{ky) - cos{nkz))]roo, 


(17) 




Hb{k) = 


/i6(e*^^'’+e /26(e*®2*’+e *®2<>e *^‘') /ate \ 


hb + 


hb + /aftC 

0 


—ikx 


V + /3(,e“*^’' /26(e*®^'’ + e-»^'2!,g*fex^ /ib(e*®'^'’ + 

-1 + 8(3 + cos{kx) - cos{ky) - cos{nkz))]Too, 




(18) 


and 




ffab(k) = 


/i(e*®i+e-*^ie**«') /2(e*®2 + _|_/'gife» \ 


3 — 201 


/i(' 

/2(e“*^" +e*^"e**-) 


g20ig iky ^ 


0 


/a + /se 


’/ g^ikx 


/3 + /^e-*"» 
0 


a — 

/i(e-'''^ 


/2(e *^2_^g*02g 


+ e*®ie*'=«) 


^ ky) 

+ (1 + i sin(n/c 


+ 


/ ^ — ikni 


/2(e*^^+e *^2g«/Ca;^ 


O i^\ fD iky 


/4COs(fc^))roo, 


(19) 


where Too = / 4 x 4 is a 4x4 identity matrix, fia, / 2 a, /sa, 
hb, hb, hb, hb, h, h, h, h, and are real coef- 
ficients, 6 >ia, 6 ^ 20 , ^ 16 , ^ 26 , and 6^2 are real parameters 
which break the TR symmetry, and n is a positive integer. 
We further define A = fs-\- and B = — f^. Different 

topological phases can be obtained by tuning the values 
of A and B in the above model. Surface energy dispersion 
can be calculated for a slab configuration and is shown 
along the momentum lines X-T and X-M in Fig. [^with a 
fixed H = 1.7 and different values of B. When B = 0.2, 
the system is a trivial insulator [Fig.j^a)]. With increas¬ 
ing H, the band gap closes on F-X between two bands 
with glide parity g^_{kx) (when B « 0.27) and drive the 
system into a phase with yi = — 1 and X 2 = 0. Energy 
dispersion for a slab configuration is shown along X-F 
and X-M in Fig. [^b) with B = 0.4. This is a semi¬ 
metal phase with a pair of Weyl nodes, emerging on F-X 
and moving in opposite direction perpendicular to f-X, 
as schematically shown by red points in Fig. if)- The 
band structure along the line connecting these two Weyl 
nodes is shown in Fig. [^g) (with B = 0.4), calculated 
by iterative Green function method. Figure [^h) shows 
the density of states in the momentum space with the 
energy determined by the nodal points, and a Fermi arc 


connects these two Weyl nodes. By further increasing 
H, the band gap closes (when B ^ 0.45) and reopens 
on X-M. One Dirac cone appears at X on the surface 
[see Fig. |^c) with B = 0.5]. This is exactly the TCI 
phase with HGC yi = — 1 and y 2 = 1. Thus, the Weyl 
semi-metal phase separates the TCI phase and the trivial 
insulating phase in our model. Another topological phase 
transition occurs for an even larger B and the system first 
evolves to another Weyl semi-metal phase with xi = ~2 
and y 2 = 1 (when 0.55 < B < 0.64) and then to a TCI 
phase with xi = and y 2 = 2 (when 0.64 < H < 1.4). 
The corresponding surface bands from the calculations of 
a slab configuration are shown in Figs, [^d) and (e) for 
the Weyl semi-metal phase {B = 0.59) and TCI phase 
{B = 0.8), respectively. 

To check if these gapless Dirac surface states are stable 
or not, we add an additional 2D layer, which preserves 
the symmetry group pgg, on the top surface of the slab, 
and introduce the coupling between this 2D layer and the 
top layer of the slab. By tuning the parameters for the 
model of this additional 2D layer, we can move energy 
bands with the required Irreps in this 2 D layer to the bulk 
band gap and couple them to topological surface bands. 
We explore the case with two Dirac cones {B = 0.9), cor- 
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r ^— X —^ M 



FIG. 9. (Color online) Band dispersions around X of the 
model of TCI in the pgg group in a slab configuration. 1.7 
is fixed and B is tuned from 0.2 to 0.8 in (a)-(e). We only 
show the surface states on one surface of the slab. The red 
(blue) lines denote the bands with even (odd) glide parity, (f) 
Schematic plot of the positions of Weyl nodes (marked by the 
red dots) in the surface BZ when B = 0.4. (g) The density 
of states calculated by iterative Green functions on the line 
connecting two Weyl nodes with B = 0.4 in a semi-infinite 
conhguration. (h) The density of states of the Fermi surface 
passing through the nodal points with B = 0.4. The other 
parameters in all of the figures are /la = / 2 a = — fib = f 2 b = 
0.1, ha = fL = hb = fL = 0 , Ola = 02a = ^16 = ^26 = 0, 
h = 1, U = 0.1, 0i = 7 r/ 2 , 6>2 = 0 and n = 2. 





FIG. 10. (Color online) The band structure of the slab with 
a 2D layer on the top surface on Y-F-X. The original slab is 
in xi — “2 phase, (a) All the bands of the 2 D layer belong 

to Irreps Ai at f. (b) All the bands of the 2D layer belong 

to Irreps Bi at f. (c) The bands of the 2D layer belonging 

to Irreps Ai and Bi coexist at f. (d) The surface states are 

gapped by adiabatically tuning the parameters in (c). 


responding to the TCI phase with xi = —2 on X-F-Y. 
When the bands of the 2D layer belong to Irreps Ai (or 
Bi) at F, the gapless surface states remains stable, as 
shown in Figs. [Tol^a) [or (b)]. However, if the bands in 
Irreps Ai and Bi coexist at F, we find a local stability 
for gapless surface bands, as shown in Fig. [T^c). The 
crossing on F-X between the bands in Irreps Ai and B 2 
at F are protected locally by glide plane symmetry. By 
tuning the parameters adiabatically, this crossing point 
can move from F-X to Y-F, leading to an anti-crossing 
between these two bands due to the same glide parity, as 
shown in Fig. 10 [d). This leads to a trivial phase. Thus, 
in this case, only one surface Dirac cone is stable and 
the topological classification is Z 2 . Therefore, our calcu¬ 
lations based on an explicit tight-binding model indeed 
support our classification of TCIs in the pgg group. 

In the spinful case, the representation of g^^y) 
is ® Di/ 2 {gx(y)) = 0 

Di/ 2 {'mx{y)) = , where D 1/2 is the 

transformation matrix of a spin-| spinor and the non¬ 
primitive translation has no effect on spin. As discussed 
in details in the Appendix spin degree of freedom in¬ 
troduces an additional coefficient for the factor system of 
the projective representations. The representations sat¬ 
isfy 


D{mh)' D{myy = asD{myy D{mh)' (20) 

where as = <^<^ 1 / 2 - Since rrixTriy = C 2 (z) and myirix = 
Q{y)C 2 {z), where C 2 {z) denotes the tt rotation around 
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the 2 ; axis and Q{y) is for the 27r rotation around the y 
axis, 0 ^ 1/2 = —1. Then as = aa^j^ = 1 at X and Y, 
and as = — 1 at T and M. Thus, all the states at T 
and M, instead of X and Y, are doubly degenerate. The 
topological classification is similar, but now depends on 
the Irreps at X and Y. 

Now, let us discuss the possible TCIs in the 'pmg 
group. There is one mirror operation and one glide op¬ 
eration in the pmg group, for which we can take rrix and 
9y — {^y\^x} with Tx = (1/2,0) as an example. The fac¬ 
tor group of the pmg group is the point group C 2 v^ which 
is generated by two mirror reflection symmetries with 
perpendicular mirror axes. The high-symmetry points 
X = (tt, 0) and M = (tt, tt) in the surface BZ possess 2D 
Irreps in the spinless case, due to the anti-commutation 
relation between mx and gy at these two points, as illus¬ 
trated in Ref. [69]. These two points are connected by 
a MIL. The coexistence of two double-degenerate points 
and the MIL which connects these two points leads to a 
mixed type-I-II TCI phase. The topological property of 
the pmg group can be characterized by the HMC defined 
on the HMP which projects onto the X-M line and the 
topological classification is Z, more than the Z 2 classifi¬ 
cation discussed in Ref. [69|. The points X and M can 
also be connected by the HSL X-T-Y-M, where T-Y is 
invariant under X-T and Y-M are invariant under 
gy. The states at T and Y belong to ID Irreps of the 
group C 2 V Similar to the case of the pgg group, if all 
the states at T and Y belong to Ai (or Bi), the states 
in the line X-T-Y-M can be labeled by the parity on one 
segment, such as X-T. This allows us to define an inte¬ 
ger topological invariant with the same form as Eq. (12), 
which is a mixture of HMC and HGC and leads to a Z 
classification. Therefore, depending on the Irreps at T 
and Y, we can have either Z^ or Z classification for the 
spinless case of the pmg group. 

In the spinful case, similar to the analysis of the pgg 
group, the spin degree of freedom gives 0 ^ 1/2 = —1 on the 
four high-symmetry points. Thus, as = <a<ai /2 = I at X 
and M, and = — I at T and Y. In this case, the double 
degeneracy occurs at T and Y. Nevertheless, topological 
classification is still similar as that in the spinless case. 

In the pAg group, there are a (^4 rotation symmetry 
and four glide symmetries gx = gy = {niylr}, 

Qd = {md\T} and g'^ = where r = ( 1 / 2 , 1 / 2 ), 

md transforms (x^y) to (y^x) and transforms (x^y) 
to {—y, —x). In the surface BZ (see Fig. CTb)), the point 
group is C^y at T = (0, 0) and M = (tt, tt), and is C 2 V at 
X = ( 7 r, 0 ) (Y = ( 0 , 7 r) is equivalent to X in this group). 
There are three independent glide invariant lines T-X, X- 
M, and M-T. The projective representations of mx and 
my satisfy 


D{mx)D{my) = a{k)D{my)D{mx), ( 21 ) 

where a(k) = e^(k-m.k).Xg-i(k-m,k).x^ X point, 

a = —1, while at T and M, a = 1. The glide plane sym¬ 
metries guarantee the double degeneracy at the X point. 


It can also be checked that for any two symmetry oper¬ 
ations in the C^y group, denoted as a and 6 , if [a, 6 ] = 0 , 
we have a = = I at T and M for the pAg group. 

Thus, the representations at these two points are conven¬ 
tional representations of the C 4 y group. A more rigorous 
proof is shown in Appendix The C^y group has both 
the ID and 2 D Irreps. If the surface bands belong to the 
ID Irreps at both T and M, no TCIs can exist. If the 
surface bands belong to the ID Irrep at T, but belong to 
the 2 D Irreps at M, the HGC can be defined along the 
line X-M, yielding the mixed type-I-H TCI phase. An¬ 
other independent HGC is possible to exist along the line 
X-T-M, depending on the ID Irreps of the states at the 
momentum T. Here, X-T is invariant under gy and T-M 
invariant under gd- At T, the glide symmetries gy and gd 
behave the same as my and md^ which satisfy the rela¬ 
tion md = C^my, where C 4 is the fourfold rotation in the 
anticlockwise direction and transforms {x,y) to {—y,x). 
The group C 4 y has four ID Irreps and Bi {i = 1,2). 
In the Irreps Ai, the character of C 4 is I, and the mirror 
parities of my and md are the same, while in the Irreps 
Bi, the character of C 4 is — I, and the mirror parities of 
my and md are opposite. This is similar to the case of 
the pgg group. If all the bands at T belong to the Irreps 
Ai (or Bi), the HGC can be defined on X-T-M, leading to 
the Z X Z classification. Similar analysis can be applied 
to the case with 2D Irrep at M and ID Irreps at T. If 
the surface bands belong to the 2 D Irreps at both T and 
M, three independent HGCs can be defined on the lines 
T-X, T-M, and X-M, giving rise to a Z^ classification. A 
model of the pAg group in the spinless case is given in 
Appendix 

In the spinful case, spin gives an additional coeffi¬ 
cient ai /2 = —I at T, X, and M. For T and M, 
as = aaii 2 = — 1 , the spinor representations belong to 
class Ki and all the Irreps are 2D. Thus, all the bands are 
doubly degenerate at T and M. For X, we have as = 1, 
so the spinor representations at X belong to the class Kq 
(conventional representations of group C 2 y)- All the con¬ 
ventional Irreps of C 2 V are ID, and thus no degeneracy 
occurs at X. One HGC can be defined on the glide invari¬ 
ant line T-M. Another possible HGC on the line T-X-M 
depends on which Irreps the bands belong to at X, which 
is similar to the case of the pgg group. If all the states at 
X belong to the Irreps Ai (or Bi) (i = I, 2) of group C 2 v^ 
the HGC can be defined on the line T-X-M, and leads to 
a Z X Z classification. 


IV. CONCLUSION 

In this paper, we have developed a theory to system¬ 
atically classify TCI phases based on the representation 
theory of 2 D space groups. We have shown that the clas¬ 
sification sensitively depends on the Irreps of the states 
at certain HSPs and HSLs. Our theory provides a ba¬ 
sis for the search of realistic materials for TCI phases. 
Since our theory is based on a semi-infinite system with 
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a specific surface, one can first identify which type of 
3D crystals can allow for the surface with the required 
2D space group symmetry. This can be achieved with 
the help of the discussion of the layer groups and the 
scanning tables in the International Tables for Crystal¬ 
lography m, as discussed in details in Appendix B of 
Ref. [69]. For example, the diamond and spinel struc¬ 
tures are described by the 3D space group Fd3m (227), 
for which the ( 110 ) surface possesses the symmetry of 
the pmg group. Therefore, one can look at the specific 
surfaces of the corresponding materials with the help of 
some crystal databases mi Eg. In addition, topologi¬ 
cal phases have also been discussed recently in cold atom 
systems with optical lattices mHzz] and photonic crys¬ 
tal systems [78H83] . Our classification, as well as the toy 
model, can also help to design crystal structures with 
the required symmetry in these systems. Disorder ef¬ 
fect can break the crystalline symmetry of the system 
[39l|84]. However, as discussed in Refs. [29l 1841486] . crys¬ 
talline topological phases will remain stable as long as the 
corresponding symmetry protection is preserved on aver¬ 
age. Aside from the crystalline symmetry, one may also 
ask if time-reversal symmetry is compatible to the TCI 
phases discussed above and if the combination of time- 
reversal symmetry and crystalline symmetry can lead to 
new TCI phases. For the first question, the existence of 
TR symmetry may or may not change topological clas¬ 
sification of TCIs. Take the pm group as an example. 
In the spinless case, for any state \ 2 p) on the MIPs, the 
state ©IV^) has the same mirror parity with and thus 
the MCNs will always be zero, and no topologically non¬ 
trivial phase can exist. While, in the spinful case, &\ 2 p) 
has the opposite mirror parity with \ 2 p)^ and as a con¬ 
sequence, MCNs can be non-zero and TR symmetry is 
compatible to TCI phases. The combination of time re¬ 
versal symmetry and crystalline symmetry can also lead 
to the so-called magnetic crystalline symmetry group, in 
which new topological phases indeed exist, as discussed 
in Ref. [56]. The classification of TCIs with time-reversal 
symmetry is beyond the scope of this paper and will be 
left to future work. It should be noted that the classi¬ 
fication approach in this paper is only suitable for TCI 
phases with nontrivial surface states. The classification 
based on the 2D space group is a subgroup of the 3D 
bulk topology. For a given 3D bulk, the surfaces along 
different directions will possess different 2D space groups 
with different topological class according to our results 
and the classification of 3D bulk topology in principle 
should include all possible topological invariants for dif¬ 
ferent surfaces (of course many of them give the same 
bulk topological invariant). Our classification may miss 
some topological invariants since some symmetry cannot 
be preserved by any surface (e.g., inversion symmetric 
topological insulators IHU). From the practical view, our 
classification is useful because non-trivial surface states 
are the main physical consequence of bulk topology. 
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Appendix A: A review of representation theory of 
symmetry groups 

In this section, we will give a short review of the rep¬ 
resentation theory of symmetry groups in the description 
of electronic band structures. 

Due to the periodic lattice structure, electronic states 
in a crystal form energy bands and are labeled by the 
crystal momentum k, which form the BZ. From the view 
of group theory, the crystal momentum k also labels the 
Irrep of the translation subgroup of the space group. At 
each k, we can define the wave vector group or little 
group Gk, which contains all the elements of the space 
group that leave k unchanged or map it onto an equiva¬ 
lent vector k + G, where G is a reciprocal lattice vector. 
Electronic states at each k can be described by the Ir- 
reps of the corresponding wave-vector group. In particu¬ 
lar, the dimension of the Irreps of the wave-vector group 
determines the degeneracy of electronic states at the cor¬ 
responding momentum. The representations of a wave- 
vector group can be constructed from the representations 
of the corresponding point-group that contains all the 
point group operations of the wave-vector group. Eor a 
space-group operation g = {r|R + T}, the representation 
of the wave-vector group ar k is D^{g) = 
where D{r) is the representation of the corresponding 
point group. 

Let gi = + r^}, we have 

D(ri)D{r2) = D^(gi)D^{g2) 

^ gik.(Ri+Tl)gik.(R2+T2)^k(-^^^^) 

^ gik-(R 2 -riR 2 )gik-(x 2 -riX 2 )L)(rir 2 ) 

_ ehk-r-^k).R2gi(k-r-^k).X2^(^^^^)^ 

Since k —r^^k = G, k)-R 2 _ i However, if G 7 ^ 

0 and T 2 7 ^ 0 , the coefficient co’(ri, r 2 ) = ^k) T 2 ^ 2 , 

and D(ri)D(r 2 ) = Cc;(ri, r 2 )D(rir 2 ). The existence of the 
factor Ci;(ri,r 2 ) shows that D{r) is may not be the con¬ 
ventional representations of the point group which always 
satisfies D{ri)D{r 2 ) = D(rir 2 ). Thus, D{r) is usually 
dubbed a projective representation of the point group 
belonging to the factor system oj{ri^r 2 ) [58]. It should 
be noted that oj{ri^r 2 ) depends on the momentum k. 
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The factor system is specified by h‘^ coefficients r2), 
where h is the order of the point group. The factor sys¬ 
tems can be classified into different classes. If D{r) is a 
projective representation of the factor system co’(ri,r2), 
the representation D'{r) = D{r)/u{r), where u{r) is 
an arbitrary single-valued function and |r^(r)| = 1, is a 
projective representation belonging to the factor system 
w'(ri,r 2 ) = system oj'(ri,r 2 ) 

is said to be projectively equivalent with co’(ri,r2). The 
set of all projectively equivalent factor systems is called 
a class of factor systems. 

For every pair of commuting element a and b in the 
point group, if , the two factor systems cj' 

and uj belong to the same class; otherwise, they belong 
to different classes. There is a class called Kq for every 
group, which contains a factor system with all the coef¬ 
ficients Ci;(ri,r2) = 1. The representations belonging to 
the class Kq are the conventional representations of the 
point group. Any factor system with uj{a^b) = uj(b^a) 
for commuting a and b also belongs to class Kq. It can 
be proved that for a point group that has more than one 
classes of factor systems, the ID Irreps can only exist 
in class Kq, that is, there are no ID Irreps in any other 
classes Kp ^ Kq. 

When spin of the electrons is taken into account, the 
wave functions have the form where I and 

2 denote the two spin states of the electrons with the 2 ; 
projection The spinors with spin-| transform ac¬ 
cording to Vii 2 {ro) = where a = {(Tx^cFy^crz) 

are the Pauli matrices, ft is the primitive vector along 
the rotation axis, and 0 is the angle of rotation around 
n. It has the property that Di/ 2 (^ 0 + 27 r) = ~^i/ 2 (^ 0 )- 
For two point-group operations ri and r2, the produc¬ 
tion rir2 can be expressed as a rotation around an axis 
with angle 0 (maybe with an inversion). Then the spinor 
representations satisfy 

-Di/2(ri)-Di/2(r2) = Wi/2(rir2)Di/2(rir2), (A2) 

where the factor C(;i/ 2 (^i,'^ 2 ) = I if 0 < 0 < 27r, or —I if 
27r < 0 < An. Thus, the spinor representations of point 
groups can be viewed as projective representations. 

Under the space-group symmetry operation, both the 
spatial and spin parts of the wave functions transform as 

^ ^ ^ [-^1/2 (^)]nm'0ji‘n 

j n(=l, 2 ) 

~ ^ \jn^im' 4 ^jni 

where D{gy = D{g) 0 Di/ 2 {g) is the direct product of the 
representations of the spatial and spin parts. The direct 
product is usually reducible and can be decomposed into 
several Irreps. 


Combining the spatial and spin parts, we have 

D(ri)'D(r2)' = [D{ri) ^ Dy2{ri)][D{r2) ^ Dy2{r2)] 

= [D{ri)D{r 2 )] 0 [-Di/2(n)-Di/2(r2)] 

= [w(ri,r2)L>(rir2)] 0 [wi/2(n,?’2)-Di/2(rir2)] 

= Lo{ri,r 2 )LOi/ 2 {ri, r 2 )D{rir 2 )' 

= Ws(ri,r2)-D(rir2)', (A4) 

where cjs(ri, r 2 ) = ^^(ri, r 2 )cJi/ 2 (^i, ^ 2 )- Thus, the spinor 
representation D'{r) is a projective representation of the 
point group with the factor system determined by ujs- 

Here, we focus on the point groups C 2 v^ C 4 y 

and Cey, which are relevant to the 2D space group. 
There are two classes Kq and Ki for C2v^ C^y and C^y, 
while only one class Kq for Csy. For the groups C2v^ 
C/^y and C^y, all the Irreps in the class Ki are 2D. If 
a = uj{a,b)/uj{b,a) = I for every commuting pair of a 
and b ([a, b] = 0 ), the representations belong to the class 
Kq. There is a systematic way to determine which class 
the representations belong to. It depends on the specific 
properties of the group. The complete discussions are 
presented in Ref. [58] . 

We take the pgg group as an example. The factor 
group of the pgg group is C 2 V at four momenta F, M, 
X, and Y. At these momenta, we need to use the pro¬ 
jective representations to describe the states. The gen¬ 
erator of C 2 V group can chosen as rrix and rriy with 
the relations: = e, = e, and rrixTUy = myirix, 

where e is the identity operation. In this case, the rela¬ 
tions between the representations of the little group Gk 
and the projective representations of the point group are 
D{mx) and D^{gy) = D{my). 

Thus, 


^gik. 2 x,e-*k- 2 r„j^ (A5) 


where I is an identity matrix, and similarly D{my)^ = I. 
The class of the representations i determined hy a = 
^ ~ ^ representations belong 

to class Kq {Ki). 

Another interesting group is the non-symmorphic 
group is pAg, of which the factor group is C^y at f and 
M. The generator of C^y can be chosen as the four fold 
rotation C 4 and the mirror symmetry rUx with the rela¬ 
tions: C 4 = e, = e and rUxC^ = c^rrix- The mir¬ 
ror symmetry rrix corresponds to the glide symmetry 
Qx = {mx\r}. Since gl = {e|ma,T + r}, g^oCi = {md\T} 
and = {md\c^T}, the relations between the projec- 
tive representation and the representations of the space 
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group are 

D{ci)^ = I, 

D{m^f = = e^'--^-D\gl) 

_ gik-2Xg-ik-(ma;X+x) _ ^ikx 

D{m,)D{ci) = D\g,)D\ci) 

= e*k-e-*k-Z)(md) = D{md), 
DiafDim,) = D''(c3)D'=(5,)e*" 

= e^’^yDinid) = e*''»D(TO^)D(c4).(A6) 

Assuming u(c 4 ) = and u(rnx) = the 

representations D{c^y = D{c^)/u{c^) and D^rUx)' = 
D{mx)/u{mx) are projectively equivalent to D{c 4 ) and 
D(mx) and satisfy 

[D{c^)'t = al, 

[D{mx)f = I, 

DimxYDic^y = [D{c^)f Dirux)', (A7) 

where a = . At the point T and M, a = 1. Thus, 

we have proved that in the group, the representations 
at r and M belong to class Kq. 

In the spinful case, we need to further consider the spin 
part, and thus as = olol\I 2 ^ where a is for the spatial part, 
and (ai /2 = ^)/^i/ 2(^5 is for the spin part. 

Appendix B: Mirror Chern number in the pmm 
group 

In this section, we will consider the direct calculation 
of the mirror Chern number in the pmm group and show 
why mirror Chern number must be zero for the spinless 
case but it can be non-zero for the spinful case. 

As shown in the main text, the Chern number of the 
mirror plane in the mirror even (odd) subspace is defined 
as the integral of the Berry curvature of the even (odd) 
subspace in the MIP. Let us take MIP of = 0, which 
is invariant under m^, as an example. On this plane, the 
mirror Chern number in each subspace is given by 

^e(o) X [ (^x5 ^ 2 ;) 5 

JMIP 

where kz) is the y component of the Berry cur¬ 

vature of the occupied energy bands in the even (odd) 
subspace 

^e{o)ikx •) kz^ ^kz^x,e(o){^kx-) kz^ dkx^z,e(o)(,kx^ kz^ 

(B2) 

with Ax(z),e(o) is the x (z) component of Abelian Berry 
connection of of the occupied bands in even (odd) sub¬ 
space. 

Next, let us show how the mirror operator mx acts 
on FQ(^o){kx^kz). For the spinless case, since mx com¬ 
mutates with my and does not change the mirror parity 


of the eigenfunctions under m^, two subspaces are sepa¬ 
rate under the mx operation. Thus, we can discuss two 
subspaces independently. For one subspace, we can in¬ 
troduce the sewing matrix, which is defined as 

Bi^apikxikz) — ('0i( 5 ^2;5 ^) (^x5 ^2;5/^))(F 3 ) 

where i = e, o for the even or odd subspace of 
^ 7/5 y’ikx^kz^a) is the a-th eigenfunction of H{kx^kz). 
Equivalently, we have 

|' 0 i( kx-, kz-, ay = ^ Bj^a^jkx^ kz)mx\yi{kx^ kz^ P)) 

/3 

(B4) 

The Abelian Berry connection can be obtained 
by taking the trace of the non-Abelian Berry 
connection, which is defined as {kx^kz) = 
—i{y{kx, kz, (^)\dki \y{kx, kz,P)). The x component of the 
non-Abelian Berry connection has the relation 

ax'^i-kxXz) = -i{'>Pi-kx,kx,a)\d-kJY^{-kx,kz,l3)) 
^oioikx-, {kx-, kz)Bp^(^hx^ k^) 

-iBa0{kx,kz)dk^Bp0. (B5) 

Then the x component of the Abelian Berry connection 
has the relation 

Ax{ kx, kz^ — Tr[( 23 ;,( kx, kz)\ 

= -Ax{kx,kz) - iTr[B{kx, kz)dk,B^kx, kz)], (B 6 ) 

where the trace is taken for all the occupied bands in the 
considered subspace. Similarly, the z component of the 
Abelian Berry connection is 

Az{-kx, kz) = Tr[az{-kx, kz)] 

= Az{kx, kz) + iTr[B{kx, kz)dk,B\kx, kz)]. (B7) 

Thus, the Berry curvature in the considered subspace has 
the property that 

F{-kx,kz) 

— ^kz^xi kx, kz) d—k^Az{ kx, kz) 
\^kz^x{kxi kz) ^kx^zikxi kz)\ 

= -F{kx,kz). (B 8 ) 

The Berry curvature in each subspace is an odd function 
of kx- Therefore, the Chern number, as an integral of 
the Berry curvature over the whole momentum space, in 
each subspace is zero. 

In the spinful case, the mirror symmetry mx changes 
the mirror parity of the eigenfunctions under my. There¬ 
fore, two subspaces of my are related to each other by 
mx in this case. As a consequence, one can follow the 
derivation above and obtain the result 

Fe{-kx, kz) = -Fo{kx, kz). (B9) 

Therefore, the Chern number in each subspace can be 
nonzero with the relation Ce = —Co and the MCN is 
Cm = ^{Ce — Co) = Ce- We conclude that the MCN 



21 


for a system with the pmm group must be zero for the 
spinless case but can be non-zero for the spinful case. 

For the non-zero MCN in the spinful case of the pmm 
group, we can further relate it to the HMC defined in 
the main text. This is because the Berry connection in 
the region kx G (—tt, 0 ) can be determined by that in the 
region kx G ( 0 , 7 r). Direct calculation gives 


Cm = 

Ce 

1 

27r 

/ dkxdky 

JMIP 

Feikxi 

kz) 




J — TT 

dk4 

j dkxF^iJvxi 
'0 

kz) + 

f dk 

J —TT 

xFeikxf 

kz)) 

-i. 

/ —TT 

dkzi 

/ dkxFQi^kxf 

'o 

kz) F 

jAk, 

,Fe{-k 

r; kz)) 


J —TT 

dk4^ 

/ dkxF^i^kx^ 
'o 

kz)- 

fAk, 

; Fq ( kx f 

kz)) 


/ —TT 

dkz J 

dkx i^F^i^kx 

0 

^kz) - 

Foikx-) 

kz)) 



— Xei^X’) ^z) Xoi^X'! ^z) (BIO) 

where Xe(o) = ^ fo ^kxFe(o){kx, h) is the inte- 

gration of the Berry curvature on half of the mirror plane 
between two lines which project onto the high symme¬ 
try points in the surface BZ in even (odd) mirror parity 
subspace. The definition of Xe{o) is actually the spinful 
version of the halved mirror chirality defined by Alexan- 
dradinata et al.[ 66 ] in the spinless case. 


Appendix C: Models of pAm, p31m and p6m in the 
spinful case 

1. pAm group in the spinful case 

Here, we construct a model of a TCI in pAm group in 
the spinful case. We consider a layered structure stacked 
along z direction, of which each layer is a square lattice 
in xy plane. There is one atom in each primitive cell. 
We consider the s orbital (instead of the p orbitals in 
the spinless model) of each atom in the simplest case. 
Assume there are two 5 orbitals on each atom, denoted 
by Si {i = A, B). The basis can be chosen as Tq = (|sa, t 
), where 

= (Cl) 

Ft 



M-—r—-X —r—-X 



M"— r —- K M -— r —- K 


FIG. 11. (Color online) Band dispersions of the model of pAm, 
p31m, and p6m group in the spinful case in a slab configura¬ 
tion. We only show the surface states on one surface of the 
slab. The red (blue) lines denote the bands with even (odd) 
mirror parity, (a) A trivial phase of the model of pAm group, 
with fAB = 0.8z; (b) a TCI phase of the model of p4?Ti group 
with fAB = 1.5z; (c) a trivial phase of the model of p31m 
group with I 2 = 0.7; (d) a TCI phase of the model of p31m 
group with I 2 = 2.5; (e) a trivial phase of the model of p6m 
group with fABi — (f) a TCI phase of the model 

of p6m group with /abi = 

where j = A, H, AB. In the basis Tq, the matrix form of 
the mirror-symmetry operation my is D{my) = ao C) i(Jy 
and that of the fourfold-rotation symmetry operation C 4 
is D{c/^) = cTo ( 8 ) If can be checked that the 

Hamiltonian satisfies 

H{k) = D{g)^H{gk)D{g), (C4) 


with i = A^B and a =t, 

The Hamiltonian has the form 


H{k) 


HA{k) HAB{k) \ 
HABiky HB{k) ) 


(C 2 ) 


where 


Hab = {tjo + tji[cos{kx) + cos{ky)] + tj2 cos{kz) 

+ tjs sin(/C;2)}cro + ifj[sm{kx)o-y - sm{ky)ax], (C3) 


where g = C 4 , my. Note that in the spinful case, although 
we only consider the s orbital, the bands are doubly de¬ 
generate at f, X, and M, and the topological classifi¬ 
cation is Z^. Set tAo = —^Bo = 23, = —tBi = — 8 , 

= —tB2 = —8, tA2, = tss = /a = /s = 0, tABO = 0-2, 
tABi = tAB 2 = 0.1, tAB 3 = lb Tuning /ab from 0 to 
1.5i, the system goes from a trivial insulating phase to 
a TCI phase with the HMCs x(r-X) = x(r-M) = 1 and 
x(X-M) = 0. The band dispersions in a slab configura¬ 
tion with /ab = 0.8i and 1.5i are shown in Fig. [T^a) 
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and (b), respectively. 


2. pSlm group in the spinful case 


Reading the character table of double group, we 
know that the basis | t) and | \) corresponds to the 
2D Irrep in the spinful case. Thus, we consider a trian¬ 
gular lattice with s orbital on each atom to build the 
model of the p31m group. The bases are chosen as 
^'o = (|syi,t), |sy44), |sb 4))^, where A and B 

denote two different s orbitals on each atom. The Hamil¬ 
tonian is 


where 


H(k) 


IlA(k) IlAB(k) \ 

ffB(k) J 


(C5) 


Hj = {tjo + tji[cos(ki) + cos(k2) + cos(fc2 - ki)] 

+ tj2 cos{kz) + ij3[sin(fci) - sin(fc2) + sin(A:2 - fci)] 

+ tj4, sin(fc 2 )}( 7 o + cos {/ji[sin(fc2) 

+ sin(fc2 - fci)] + fj2[- cos(A: 2) + cos(A:2 - ki)\}ax 
- e^^/^{fji[sm{ki) + (sin(fc2) - sin(fc2 - fci))sin 

+ /i2[cos(A:i) - (cos(A: 2) + cos(A:2 -ki)) sin (^)]}o'y 

(C6) 

where j = and ki^k 2 denote the momentum 

along the basis of the BZ bi = (\/3, — l) 7 r and b 2 = 
(0, l)27r, respectively. In the basis Tq, the matrix form 
of the mirror symmetry operation rriy is D{my) = (j^^icFy 
and that of the three-fold rotation symmetry operation C 3 
is D{cs) = ctq ( 8 ) ^(cTo By direct calculation, we 

get [m^,C 3 ] 7 ^ 0 , which indicates the double degeneracy 
at T, K and K'. Set tAo = —tBo = 23, Iai = —tsi = — 8 , 

tA2 = —tB2 = —8, tAS = tB3 = 0, fA\ = fA2 = fBl = 
fB2 = 0, tABO = 2, IabI = tAB2 = = 0, ^ABA = h 

and define h = (/abi + /AB 2 )e*"'/®/ 2 , I 2 = (/abi - 
fl’^- Fix = 0 and increase ?2 continuously. 
When I 2 = 0, it is the trivial insulating phase. When 
I 2 ^ 0.96, gap closes on T-K, and drive the system into 
to a semimetal phase. When I 2 ^ 1.9, the gap closes on T- 
K', and after that, the system becomes insulating again, 
with one surface Dirac cone at T. The band dispersions 
in a slab configuration with I 2 = 0.7 and 2.5 are shown 
in Figs, prj^c) and (d), respectively. 


3. p6m group in the spinful case 

The Irreps of Cqv (for T) and C 2 V (for M) group 
are all 2D in the spinful case, thus we still can con¬ 
sider only the s orbital of the atoms on a triangu¬ 
lar lattice. Meanwhile, from the analysis of the Csv 
group, we know that with s orbital the bands are dou¬ 
bly degenerate at K, too. The basis can be chosen as 



FIG. 12. (Color online) Band dispersions around X of the 
model of p4g group in the spinless case in a slab configura¬ 
tion. The parameter h = 0.3,1.2,0.8 in (a), (b) and (c), 
respectively. We only show the surface states on one surface 
of the slab. The red (blue) lines denote the bands with even 
(odd) glide parity, (d) The density of states calculated by it¬ 
erative Green functions on the line connecting two Weyl nodes 
(N1-N2) with h = 0.8 in a semi-infinite configuration. The 
line N1-N2 parallel to f-X with the middle point on X-M. 


^'o = (|sA,t), |sa 4), |sB,t), |sB,i))^, where A and B 
denote two different s orbitals on each atom. By set¬ 
ting the coefficients fj 2 {j = A^B^ AB) to zero in the 
Hamiltonian of p31m group, we get the Hamiltonian of 
p6m group. The matrix form of cq rotation operation 
is D{cq) = cTo G) |(\/3cro + icFz), and that of niy is still 
<Jo G) icFy. Set tAo = —tBo = 23, Iai = —tBi = ~8, 
tA2 = —tB2 = —8, tA3 = = 0, fA1 = fBl = 0, 

tABO = 1, tABi = tAB 2 = 0, Iabo = lb Tuning 
from 0 to 1, the system transforms from a triv¬ 
ial insulating phase to a TCI phase with one surface Dirac 
cone at F. The band dispersions in a slab configuration 


with Jabi = 0.5e and e are shown in Fig. 
and (f), respectively. 
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Appendix D: Model of pAg group in the spinless case 

In this section, we give a model of pAg group. An anti¬ 
ferromagnetic lattice can be constructed to realize the 
pAg group. It is consisted of two interpenetrating square 
sublattices, which are connected by a translation (|, ^, 0) 
to each other (the lattice constant in xy plane is set as 1). 
The magnetic moments in each sublattice has a ferromag¬ 
netic order perpendicular to the xy plane, while the mag¬ 
netic directions of the two sublattices are opposite. Con¬ 
sider the Px and Py orbital on each atom, the basis can be 
chosen as = (|^, 1,^^;), \A,l,Py), \A,2,px), \A,2,py), 
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|5, \B,l,py), \B,2,Px), |5,2,p^))^, where 1 and independent p orbitals. The Hamiltonian has the form 

2 denote two sublattices, and A and B denote two sets of 

H{k) = 

where 


( HA{k) HAB{k)\ 

\HabO^)^ Hb{^) ) 


Hja = {tjo + tjs cos{kz) + tj4 sin(fc2))roo + {tji cos{kx) + tj2 cos{ky)) 


(Foo + Foa) 


+{tji cos{ky) + tj2 cos{kx)) + tj5 sm{kx) sin(fcj^)r33; 


Hjb = {fjO + fj 3 COs{k^) + fj 4 Sm{kAi ^32 + {fjl COs(fcj:) + fj 2 COs(fcj;)) +^F32) 
+{fji cos{ky) + fj2 cos{kx)) ^ F 3 i+zr 32 ) sm{kx) sm{ky)Toi; 

Hjc = [aji{l + cos{k^ + ky)) + aj2{cos{kx) + cos(fcj^))] + ^^23) 

+ [aji sin{kx + ky) + aj2isin{kx) + sin(fcy))] 

+ [aj2i^ + cos{kx + ky)) +aji{cos{kx) + cos(fcy))] 

+ [aji{siT 3 .{kx) + sin(fcj^)) + aj2 sm{kx + ky) A^° 

(r-i-i - Too) 

Hjd = [bji{l + cos{kx + ky)) + bj2{cos{k:^) + cos{ky))]- - - - 

(r T r ) 

+ [bji sin(fca; + ky) + 6j2(sin(fca;) + sin(fcy))] ^ ^ —— 

-[ 6 j 2 (l + cos{kx + ky)) + bji{cos{kx) + cos(fcj^))] 

+ [bj 2 sm{kx + ky) + bji{sm{kx) + sin(fcy))] 1— 

~ + Hjb + Hjc + Hjd. {j = A, H, AB) 


(D2) 


The generator of the group can be chosen as glide- 
plane-symmetry operation py = {m^jr = (|, |)} and 
four-fold rotation C 4 . In the basis T, the matrix 
forms of the symmetry operations are D{gy) = ao ( 8 ) 
[e~^^^^‘^{cos{kx/2)ax + sm{kx/2)ay) ( 8 ) cr^], and D{c 4 ) = 
cro0[ie~'^^^^‘^{cos{kx/2)cro^sm{kx/2)az)0cry]. By direct 
calculation, it can be shown that at T, {^^, 04 } = 0 and 
at M, [gy,C 4 ] = 0. Thus, in the basis T, the bands are 
doubly degenerate at T, while not degenerate at M. Since 
we have proved that with any basis the bands are doubly 
degenerate at X, one topological invariant HGC can be 
defined on T-X. In our model, the representation of C 4 
for the bands at M depend on the value of the parame¬ 
ters of the Hamiltonian, thus in general, as discussed in 
the main text, the topological classification of our model 
is at least Z. Set tAo = ~tBo — 23, Iai = ~tBi — 85 
tA2 = —tB2 = — 8 , tAS = —tB3 = “ 8 , ^^4 = = 0, 


tABO = 1, tAB4 = Define h = (a^si + <^as 2)/2 
and I 2 = {ciABi — CiAB 2 )/{‘^i) and fix I 2 = 2li. The 
other parameters can be chosen to be zero, or some ar¬ 
bitrary numbers with small absolute values, which will 
not change the topological properties. When li < 0.5, 
the system is in a trivial insulating phase. Increase li to 
0.5, gap closes on X-M, and then drive the system into 
the semimetal phase with Weyl nodes. The band struc¬ 
ture along two of the Weyl nodes are shown in Fig.jl^d) 
with li = 0.8. When /i ^ 1, gap closes on T-X. In¬ 
creasing li further, the system goes into a topological 
nontrivial phase with x(r-X) = 1. The band dispersions 
in a slab configuration with li = 0.5,1.2, 0.8 are shown 
in Fig.[T^a), (b) and (c). 

Appendix E: Useful character tables 


We have the following character tables (Tables HI-V) 
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TABLE III. The character table of C^v single group. 



E 

2 C 3 

Srriv 

^1 

1 

1 

1 

A 2 

1 

1 

-1 

E 

2 

-1 

0 


TABLE IV. The character table of C^v single group. 


C4V 

E C 2 

2Ca 2mv 

2md 

Ai 

1 

1 

1 

1 

1 

A 2 

1 

1 

1 

-1 

-1 

Bi 

1 

1 

-1 

1 

-1 

B 2 

1 

1 

-1 

-1 

1 

E 

2 

-2 

0 

0 
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